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Al'T OPTIMAL SmGINEERING PLANNING MODEL 
FOR WATER RESOURCES DEVELOPMEL-IT 

The primairy purpose of water resources analysis is to 
provide the infomation necessary for optimal development of 
water resources systems. In recent years such analysis and 
design have been greatly helped by the use of systems analysis 
and operations research techniques. The optimal development 
of v;ater resources system consists of two basic issues: 
a) scheduling of promising projects and b) Operation policy 
of system to achieve the required objective. These two 
problems are solved using mathematical programming method. 
Scheduling problem decides which of the si3b--set of projects is 
to be selected, their sequence of construction and timing of 
introduction of projects into the system to provide the capacity. 
Operation policy of the system gives the releases from different 
projects of the syston to meet the demands specified. The 
demand for water is met with either by surface water supplies 
or ground-water supplies or both. Thus the model analyzes the 
integrated use of ground and surface supplies. 
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The locations of potential dain sites and well-field units 
are determined by the topography and geology of the basin, 
aquifer characteristics, runoff statistics and recharge potential 
of the basin. 

In the model the following assumptions are made: 

1) Time scale of introduction of a new project is chosen 
as one year. Each project selected, becomes the part of the 
system at the beginning of the year. 

2) Each project can be constructed to three different 
scales. Each scale of construction is considered as a separate 
project. 

3) Hydrology of each su!o-basiY\. is known. 

4) The present and future demands for various uses - 
irrigation, hydropower, municipal and industrial supply are 
known over a planning horizon. 

The optimization problem is formulated as a zero-one 
mixed integer programming problem. The problem is decomposed 
into a set of feasible combinations. Each feasible combination 
is solved by implicit enumeration using Little's Branch Sc Bound- 
algorithm. The optimal returns of each feasible combination is 
determined by Fullcurson ' s out-of-Kilter algorithm. Then that 
combination of projects which gives the maximim value of returns 
(present value of net benefits) is selected as the optimal sifo-set. 



The model is extended to incorporate stochastic environ- 
ment talcing inf lovers to the system as random variables. Chance- 
constrained approach is adopted for the stochastic formulation. 
Using linear decision rule of Revelle et al. and probability 
distribution of stream flovrs, optimisation problem is solved, 
for the confidence level as specified Ijy management. 

Both evaporation losses from the surface of reservoirs 
and return flov/s from irrigated areas and municipal and indus- 
trial supplies are included in deterministic and stochastic 
models. The return flow coefficients of 0.40 for irrigation 
and 0.70 for li & I supplies are adopted. The evaporation losses 
are taken as constant in the season of the year. 

The model has been proved as an efficient technique for 
scheduling and sequencing of projects and in determining effi- 
cient operating policy. The model is flexible to incorporate 
any institutional constraints lilco contingent and compound 
pro j ects. 

Integer values which have zero-one variables indicate 
V7hich projects should be constructed, and to what capacities in 
order to meet the demands. Activity levels of continuous 
variables indicate the quantities of water released or pumpages 
from the projects of the system. 

The model is formulated for the seasonal storage require- 
ments. In the model the year is divided into two sessions - 
Khar if (May to October) and Rabi (November to April) . 
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With access to a large enough computer it is possible to 
use this model over many time periods depending upon the charac- 
teristics of inputs and demands of the basin. 

The methodology is applied to the planning process of 
Cauvery river basin (Karnataka, India) , The basin under consi- 
deration drains an area of 36240 sq. km. The main physiographic 
divisions of the basin are (a) VJestern Ghats (b) Plateu of 
Karnataka, The delta region of basin which is in Tamilnadu is 
not considered in the analysis. However, the mandatory releases 
for down stream uses are incorporated by suitable constraints. 

The basin consists of a number of tributaries. There are 
24 potential dam sites in addition to the existing 12 reservoirs. 
There is a provision for using ground water for supplementing 
irrigation in certain regions of the basin. The demand comprise 
of irrigation, hydropov;er and municipal and industrial supplies. 

At each location of the project (either surface v/ater or 
groundwater) , three scales of development are considered and 
each scale is considered as a separate project- For each project 
capital cost, operation, maintenance and replacement costs and 
annual maximum returns are considered. 

The problem is solved on IBM 7044/1401 system. The results 
indicate the site and timing of projects and the operation of the 

system for maximizing present value of net benefits. 
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Sensitivity analysis v&s conducted for different parameters 
namely discount rates, time horizon, and capital costs of the 
projects. It was found that the time horizon of more than 30 
years has no effect on the scheduling and sequencing of projects- 
Social discount rates of 3 %- 5 % are used in the model and results 
are presented. 

For the successful application of the model a considerable 
input from economists, engineers and other experts familiar vjith 
local situations is required. It can be concluded that the model 
serves as a powerful methodology for the investment policy 
iDroblems. 

The model -can be extended for stochastic demand pattern on 
the system. 



Chapter 1 


INTRODUCTION 


1 . 0 General 

For a long period, conventional methods were employed 
in the analysis of water resources system. In studying river 
basins, v/hich are enormously complex, various interacting 

physical, economic, social and political systems should be 

-ter 

talcen into consideration. At the outset, systems study of a 
river basin, an, approach is requiired which will reduce this 
complexity to the level at v/hich the systans analyst can perceive 
the major relevant interactions in that basin. The appropriate- 
ness of an approach depends upon the systems analyst and his 
previous experience and the nature of physical system under 
consideration. 

In recent years there has been a growing interest on the 
part of government agencies, planners, engineers and economists 
in improving the methods of analysis of water resources systems. 
Much thought and v/riting has been devoted to the subject, but 
most of it is problem oriented. Many reports and research 
articles have been published which are based on special situations 
and technical conditions. 
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Water resources development is required to minimize the 
discrepancy that exists between natural supply of water in 
time it space - both in quantity and quality - and the demand 
for it. To bridge this gap in demand and supply it is 
necessary to construct hydraulic structures and specify their 
operational policy. 

Water resources systems include physical/ economic, 
political, legal and ecological aspects. It is necessary to 
identify and integrate all these aspects in order to achieve 
comprehensive development of v^ater resources - 

In the present study an attempt is made to use systems 
analysis to analyze alternate investment patterns in water 
resources projects and to present efficient operation policy. 
The result of this v/orlc is a nev7 planning model based on branch 
and bound algorithm combined xvith out-of-Kilter algorithm for 
schedijling and sequencing of water resources projects. Use of 
this method is demonstrated by applying it to the Cauvery river 
basin. It is assumed that hydrology of the basin and other 
institutional and "oudgetary constraints are hno^m. 

1.1 Planning Horizon 

A review of the performance of existing vrater resources 
projects shows that 10-15 years or more elapse between the 
conception of a need and the full utilization of the system. 
Apart from political hassal period in which, rights, responsi- 
bilities and obligations of all concerned are thrashed out in a 



political envircnment/ present day water resources systems 
require 5-10 years or more as construction period. Water 
resources systems are quits duraJsle. Systans that were built 
centuries ago are still in use in many parts of the countiry. 
-■'iowever, for economic analysis, a planning horizon of 30 years 
is taken. Planning horizons of 40 and 50 years are also 
examined. 

1.2 Discount P.ate 

The classical economic theory views interest rate as a 
rate of discount on future goods and services whose level is 
determined by balance of considerations of time preference and 
productivity (Mao, 1969). 

The concept of interest rate is of considerable abstrac- 
tion since it is observed that in the market a n'omber of 

different rates vary primarily with length of time and type of 
security offered. The theoretical interest rate refers to 
rislcless investraents and may vary with time span. 

In case of water resources system there is another 
important factor to be considered. Huge investments are made 
during the planning and construction period (Scl.stein, 1965). 
The returns start accruing only after the project is completed 
or has attained a certain stage of maturity. This major 
difference in time betw^een the centre of gravity of expenditur 
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and the centre gravity of returns leads to a requirement of 
a means of coraparing the value of money at different times. 

The accepted method is to calculate the discounted 
return according to som.e discount rate. VThat the rate should 
be, is still a subject of controversy. With long horizons 
involved, low discount rates aid the feasibility of project 
while high rates can easily kill it. 

When applied to v/atsr resources systems analysis, low 
discount rates are justified. In this study discount rates 
of 3-5 percent are considered. 


1,3 Hydrology 

Hydrological variations create a need for dams and 
reservoirs in order to regulate the stream, flov; to make it 
correspond more nearly to the demands im.posed upon it. This 
regulation may be for the purpose of controlling high flows 
or for auejmenting the Icrw flov^s or for both. 

Stream flow is essentially a stochastic quantity. 
Regardless of regulatory storaeje provided or manner of its 
release, there still exists a certain probability that the 
high flov. 7 s will not be contained or that there will be insuf- 
ficient water to augment the low flow. 

In the stochastic model of this study, the stochastic 
aspect is considered and the chance constraint approach is ■ 
used in the analysis. 
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Evaporation losses are quite considerable in arid and 
semi-arid regions. It is estimated that 1.2 m to 1.5 m are 
lost due to evaporation from the water surface of reservoirs. 
These are given due consideration in the analysis. 

1 . 4 Water Demands 

The demand for water is of tvro types - consumptive use 
and non-consumptive use. The consumptive use is generally for 
agriculture, municipal and industrial supplies. Non- consumptive 
uses include recreatiori, navigation, pov/er generation and wild 
life preservation. 

The demands of water are obtained by knowing the irrigable 
area, type of crops, farm technology'' adopted, population growth, 
per capita consumption and standards of living. These demands 
are estimated sub-basinv/ise in the entire basin. 

The best indicators of the demand pattern over the 
intem/al - say a year, is the maximim water demand, 

minimum demand, D . , The total demand over the interval is 

mm 

the area under the demand curve, y = d(t), figure 1. 

T_.. 

D = i d (t) . dt 

o J 

D is talcen in the units of hectare meters. 
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Tic. 1.1. Demand Ciorve 

1.5 Cost Characteristics of Water Resources System Elements 

The cost functions associated with various components 
display ‘economy of scale’ to a greater or lesser extent, i.e,, 
cost per unit of capacity decreases as the capacity increases 
(Hall Si Dracup, 1970; Hirshliefer et al., I960), There is 
nearly alv/ays a finite discontinuity of the cost functions at 
zero capacity. In addition, most components of the system 
involve very large investments compared with operation and 
maintenance costs. It is rather difficult to derive generalized 
expressions for cost functions of dams. At any site there will 
be su]3stantial fixed costs associated v/ith preparation of site, 
for construction particularly for foundations and abutments. 

The relative economy of scale clearly depends on the cost of 
site preparation and height to width ratio. As pointed by 
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Hall and Dracup (1970), depending upon the dam site and scale 
of development, economy of scale could provide a doubling of 
capacity for as little as 1.2 times the cost or might reguire 
twice or more than tv/ice the cost. In many cases V 7 ith good 
sites, the physical limitation of stream flow may limit the 
height long before the cost function reaches a convex segment. 

1.6 Water Rates 

Returns from system related to water use are even more 
difficult to treat. Water rates are fixed depending upon type 
of the use. In arid and semi-arid regions, sufficient quantity 
of water has to be supplied to provide for the evapotranspira- 
tion and to leach salts, if any, accumulated at the root zone. 
Hence, irrigation requirement is a major demand for water. 

In India there exists a wide variety of methods for 
calculation of water rates supplied to agriculture. These 
rates vary from state to state, from basin to basin and from 
project to project. The value of crop at maturity served as a 
rough guide in fixing irrigation water rates (C.V-T. and P.C., 1955) 
Irrigation Rates Committee, U. P., India (1939) reports that the 
value of water is a function of cost of the supply and increase 
in the produce from the land irrigated. A fair rate is one 
which yields adequate profit on out— lay on one hand, and offers 
sufficient inducement to the cultivator to avail the facility, 
on the other hand; and also prevents them from extravagance or 
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v/aste. The national Council of Applied Economic Research 
(NCAER, 1959) suggested that whatever be the net benefit, the 
actual water rate should in no case exceed 50 percent of extra 
benefits available to cultivators. The actual charge could 
vary betv/een 20 percent to 50 percent of net benefits. 

The system of charging for v/ater according to volume 
supplied is scientific and rational in nature. But it involves 
large initial investment on measuring devices and supervise 3:^'" 
staff. Moreover, there is the rish of damage as the people may 
tamper with the meters. It is also readily susceptible to 
malpractices on the part of the officials. 

Though the net benefit method is accepted as a sound 
basis for fixing of water rates, it is rather difficult to 
apply .it as it involves collection of data on type of crop, 
yield, prices and cost of production. Generally vxater rates 
are fixed depending on the type of crop grov/n and acreage under 
different crops. This practice indirectly implies that v/ater 
rates are fixed on amount of water used, efficiency of irriga- 
tion and conveyance from head works to the fields. NCAER 
recommended that compulsory rate of 15 percent to 30 percent 
of net benefits is to be levied as xirater rates. For this 
purpose, the area would be divided into various groups and 
net benefits accruing to cultivators are measured for different 
crops in each soil and each season. Another study in this 
connection, v/orth mentioning is: Evaluation of Major Irrigation 
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Projects - Some Case Studies - Planning Commission, India (1965). 
The details of finding the water rates from the revenue collected 
is given in Appendix. a. 

Municipal and industrial water rates are talcen from the 
schedule of ’water rates applicable to the city of Bangalore. 
Similarly, the revenue of water diverted for power generation 
is obtained by tariff schedule of Karnataka Electricity Board 
(KEB) . 

1.7 Conjunctive Use of Ground and S\irface Waters 

There is a grov/ing awareness for conjunctive utilization 
of surface and ground v/ater supplies. As the increased demand 
for vjater is not m.at xfith surface supplies, the conjunctive use 
of ground and surface supplies comes into sharper focus. This 
implies an allocation problem involving both physical and 
economic parameters and this allocation must be done in an 
efficient v/ay. The required allocation is best -dealt with in 
the form of tv/o models - a physical model and an econom.ic model. 
The former enables the physical quantities of available water 
under varying cond.itioriS to be determined and the latter allows 
the allocation of this water tc be made to various demands, 
subject to given costs and the optimum allocation is deterrained 
iteratively to r.ieet the specific objective. 

In this dissertation conjunctive utilization is consider -ed 
implicitly and the model decides when and in what quantity the 
ground water is to be utilized. 
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1 . S Return FI ows 

Not all the water diverted to agriculture is consimp- 
tively used by the crops. That part X'^!'hich is not consumptively 
used runs off the crop land as surface flow or seeps into ground 
and returns to the rivers as return flow (Maass et al., 1962). 
Some part of the vratcr which seeps into the ground becomes part 
of the -water called interflow and is essentially available as 
surface water since the streams, labos or reservoirs intercept 
it. The remainder becomes part of ground water sup]ply by the 
process of deep percolation.. Return flow coefficients are 
estimated for different hydrological units and considered in 
the analysis. 

Similarly not all the water diverted for municipal and 
industrial supply is consumptively used. Waste v/ater from 
residential areas and industrial plants after treatment is 
channeled into surface stream.s and is also knovm as return 
flows. The total return flov/ from agriculture and municipal 
and industrial supplies is taken into account. 

The various parameters considered above are used in the 
development of subsequent raodels. 


1.9 

Chapter 2 provides a revievr of current practices on 
water resources planning. 
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Conventional analyses are conducted on a partial or 
project by project basis but the selection and performance of 
any specific project at any particular time is heavily influenced 
by the composition of the rest of the projects; Marglin (19S3) 
says ’Planning of water resources development projects is more, 
even more complex by the impossibility of separating the problem 
or choosing the optimal input-output bundle from the problem of 
scheduling the construction of the component increments- We 
cannot, for example, simply determine the optimal input-outout 
bundle as if the entire plan is undertaken today and than 
determine the optimal sequence of construction of various 
phases of the plan. It is often the case that input-output 
bundle v/hich v/ould be optimal under the assumption that all 
construction m.ust take place at once is no longer optimal when 
we remove this restrictive assimption and expand our horizons 
to include future construction possibilities . . - ’ . 

In the past, systems analysis techniques could not be 
used extensively because of the time and money required for the 
analysis considering all the components of the system. Hov/ever, 
with computing facilities becoming v/idely available, this is no 
longer a problem. ’■'fnen huge investments are involved in v/ater 
resources development it is liic'hly desirable to adopt systems 
analysis and operations research techniques. 

Chapter 3 presents different m.odels v/hich are developed 
for the sequencing and scheduling water resources projects and 
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to find out efficient operation policies. Any model developed 
to solve v;ater resources problems should have the following 
characteristics: 

1. The model must optimize som.e specified objective. 

2. The model should be feasible. It should be possible to 
set up and solve the model within a reasonable amount 
of time. 

3. The m.odel should be efficient in data use. It should be 
possible to use the relevant available information and not 
insist on data which one is not likely to obtain. 

4. The m.odel should yield such results which can be understood 
by those who actually make the final decisions and can be 
implemented. 

5. It should be flexible enough to allow for sensitivity 
analysis of important parameters. Particularly attention 
should be given to those paraiTieters v/hich are uncertain 
in nature. 

The development of the investment model in the deter- 
ministic enviroriTiient followed by stochastic model is presented 
in Chapter 3. Finally evaporation losses are also included in 
the analysis. The return flov7S are talcen into account in all 
the models. 

The model is applied to a largo river basin system - 
Cauveri'- basin in South India. The basin features are given 
in Chajpter 5 . The results and analysis are presented in 
Chapter 6 along V 7 ith the conclusions drawn. 



Chapter 2 


CURRENT PRACTICES IN 

WATER RESOURCES PLANKING AND DEVELOPMENT 


2 . 1 Introduction 

It has been stated that planning is the process by v/hich 
society directs its activities to achieve goals it regards as imp- 
portant (Veisman, 1974). Hovjever, where available water and 
related land resources are scarce and the current projected 
competition among various users is great, it becomes exceedingly 
difficult for planners to (l) properly define goals and objec- 
tives, (2) formulate plans that satisfy competing demands or 
conflicting objectives, (3) gain the puJolic acceptance of 
proposed plans, and (4; finally implement and operate a system 
of facilities so as to meet society's goals (Texas Water 
Development Board - TVJDB, 1968). 

Over the past decade, engineers, economists and systous 
analysts, to the best of their ability and within the limits 
of computational capabilities, have defined meaningful advahced 
planning procedures and supporting analytical techniques. 

Largely, their use has resulted in planning strategies that 
depend upon the fast computational capabilities of the computer 
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to evaluate plans and operational criteria identified as 
attractive by the planner in his attempt to find the optimal 
strategy/-. Terras such as minimum cost, maximum return and 
maximum net benefits have emerged and have been extensively 
used as the basis for quantifying optimality. 

Investiuent planning requires a determination of the 
goals and objectives of the government, of operating policies 
to be followed in future time period and resource requirements 
of alternative investments. 

Any systems analysis model must be designed to interact 
in two directions - with other such models and with exogenous 
agencies. Models can be classified as preliminary screening, 
detailed investment analysis and operating policy models. Flow 
of information within the systems analysis will be from the 
preliminary screening to the detailed investment analysis to 
the operating policy models v/ith feedback from the latter model 
to the detailed investment analysis model. 

There is a voluminous literature on the use of systems 
analysis and operations research techniques in the field of 
v/ater resources systems analysis. Available literature can be 
grouped under tv70 broad categories, a) planning and design of 
water resources, and b) operation policies of water resources 
system. Most of the worlr. reported in the last 10-15 years 
has been devoted to the second categori^ of problem - i.e., 
finding out operating rules of existing water resources system. 
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Although operation of single reservoir is well documented, 
there are verg' few worhs in the literature on hovj to size 
and operate multi-unit-multi-purpose v/ater resources systems. 

One of the early works and most extensive work in v/ater 
resources systems v/as carried out by Maass et al., ( 1562 ). 

They introduced objectives, defined concepts, used systems 
approach and mathematical programming methods. 

In capital investment planning of water resources 
projects, there exists both a set of durable high cost 
projects and a number of demands over some specified time 
horizon. The scheduling problem considered here is: 

(i) which subset of the projects in the basin to select, 

(ii) in what sequence to construct the selected projects, 

(iii) at v/hat time to bring each of the selected projects 
into the system, and (iv) to specify the operating procedure 
in order to meet the demand at every point in the basin and in 
the planning horizon and to optimize some prospecified objective. 

A special case of scheduling problem is the sequencing 
problem. V'Jhereas, scheduling probleras involve project selec- 
tion, sequencing and timing decisions, segi-iencing problems 
involve only project sequencing and timing decisions. 

Sequencing problem arises when, for example, the projects 
in the given set must be constructed in order to satisfy the 
requirements. In this case project selection decisions are 
obviously not necessary. 



16 


Mathematical programming methods have been applied to 
capital investment problems in water resources systans analysis. 
Marglin (1963) perforrried an economic analysis of selection of 
projects and timing of project construction. Calculus of 
variation and linear programming methods were used in this 
analysis. 

Buras (197 2) discussed the application of mathcsmatical 
programming methods to the design and operation of multi- 
structure system. hayak and Arora (1970') used separable 
programming technique for finding the minimumt capacities of 
reservoirs in a system designed for flood control, low flow 
augmentation and recreation. Young, Mosely and Evenson (1969) 
used pattern search technique for sequencing of reservoirs. 

Many dynamic programming approaches vjere presented for sequencing 
of projects (Morin, 1973; Morin and Esogbue, 1971). Butcher, 
Haimes and Hall (1969) presented a dynamic programmdng algorithm 
to solve a special case of the sequencing problem for certain- 
water supply projects. Burton (1959) formulated an integer 
programming model of project relation aspect of the capital 
investment problem. Inventory theory was applied by Manne 
(1962, 1967) to project timing and project capacity decisions 
encountered in capacity expansion for developing nations. 
Mulvihill and Dracup (197 4') developed multi-level solution 
technique to determine the optimum sizes, timing of a conjunctive 
urban v/ater supply and waste v;ater system. Trott & Yeh (1973) 
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and Becker & Yeh (1974) applied dynamic programming approach 
to Eel river basin - Sacrsmanto, California - for finding 
optimal sizing and sequencing of multiple reservoirs. But 
Morin (1973) pointed out that conventional dynamic programming 
may yield non-optirnal solution to sequencing problem, Orlob 
(197 0) discussed the approach taken by planners for Texas v/ater 
system. Trans Texas Division of Texas water system vrould 
comprise of IS reservoirs, more than 800 Imi of canal networks. 

In addition, there xfould be pumping facilities to raise the 
water from sea-level to over 1000 m elevation. Their planning 
problem was: 

Given; (l) The location of all reservoirs (2) details of canal 
net-work (3) schedules for in-basin demand for each reservoir 
and each junction of the system (4') cost of construction, 
operation and maintenance of all elements (5) hydrology of supply 
for each reservoir. 

Find: Optimal schedule of projects to meet the specified 

demands within the prescribed legal, institutional and physical 
constraints . 

The Texas Water resources group used a coirilDination of 
linear programming, simulation technique, response surface 
methods. This method was to seek near optimal solutions 
rather than exact optimal solution to overcome the limits 
that existed on computation time and computer facilities. 
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Simulation techniques are applied in the analysis of 
water resources systans (Maass et al., 1962'), Simulation 
process attempts to reproduce the system as faithfully as 
resources permit. The computation permits the nonlinear, 
dynamic, stochastic and feedbaclc responses of the system given 
all design and operating policies before hand. But the cost 
of simulation in case of water resources systems involving 
extensive adjustraent of parameters mahes this technique 
prohibitive. 

Kufschmidt and Fiering (1966) carried out the systematic 
simulation study for optimal sizing of reservoirs in a water 
resources system designed to maximize the return on investment 
and to meet the schedule of v/ater demands. In the analysis 
it was assumed that the configuration of dams and operating 
policy for each dam and hydrology were known. For each combi- 
nation of dam sizes the simulation gave the revenue from 
operating the system. Simulation technique was also used by 
VJeiss and Beard (1972) to analyse the water resource problem, 
namely, (l) which of the specified set of projects should be 
constructed, (2) how large each of the projects should be, 

(3) timing of construction, (4) how should the systen be 
operated so as to minimize the capital cost plus operational 
costs over the planning period. 

Techniques of optimization and simulation v/ere merged 
in the analysis presented by Viessman et al., (1975) to select 
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the most efficient arrangement of components for regional 
water resources development and management. 

O'Laoghaire and Heramelblau (1971 ) presented method for 
optimal capital investment in expansion of an existing water 
resources system. New approach vjas presented using zero-one 
programming and networlc flow algox'it'nras . Deterministic 
hydrology x-.’^as used in the analysis for sequencing of surface 
supply projects. 

/Ground water aquifers have also been long recognized 
as important sources of water. Hov/ever/ in the past no consi- 
deration was given for the interaction heti-aeen surface and 
ground waters and to establish a rational basis for the 
development and use of subsurface storage in water resources 
development. In the last decade many investigators have given 
the operating policies for conjunctive use of surface and sub- 
surface supplies. The concept of conjunctive use of all 
sources of water has been developed in recent years (Tyson and 
Weber, 1964; Koenig/ 1363, Dracup, 1966; Rogers and Smith, 
1970). 

Chun et al., (1964) formulated alternative plans to meet 
the estimated schedule of future water demand over a planning 
horizon of 25 years. The demands could have been met entirely 
through an all-surface distribution system or by all ground 
water system or by a combination of these. 
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Dracup (1956) used linear programming method to find 
optimum policies for conjunctive use of ground and surface 
waters. Hatfield and Pattison (1970) applied mathematical 
prograrriming approach for optimal development of surface and 
ground supplies in Hunter Valley Easin^ Australia. A compre- 
nensive review on conjunctive use v/as given by Milligan (1970) . 
A chance constrained method was used by Smith (1970) for the 
conjunctive use of surface and ground v;ater supplies to meet 
the irrigation requirements over a planning horizon. 

For a given set of hydrological data and aquifer 
characteristics/ the optimization of the solution of the 
conjunctive use system involves the solution of three problems: 

(a) deterrai nation of design capacities for surface facilities 
includinc; recharge facilities, 

(b) deternination of the system service area, 

(c) determination of operating policy specifying reservoir 
releases and aquifer pumpages. 

There is a voluminous literature on the operation 
procedures of reservoirs. (Revelle et al., 19S9; Nayalc and 
•^^ora, 197 0). Most of the literature available on investment 

policies in v/ater resources systeras deals v/ith a specific 
economic situation or a particular institutional constraint. 

Additional difficulty with the literature available in 
this field is that it originates from advanced countries v^hero 
the economic and institutional constraints are quite different, 
from those of the developing countries. 
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From the brief simmary of literature review it can be 
concluded that the systems approach has been used in the field 
of water resources engineering and problems have been solved 
using programming methods and simulation techniques. However, 
only a few investigators have considered the random nature of 
inflov7s and demands of water. Horeover, evaporation losses 
v/hich are important especially in arid and semi-arid region 
in the efficient management of water resources have not been 
given due attention in the optimization m.odels. 

Here an attempt is made to formulate the optimization 
model for investment decisions and to obtain operating rules 
for different components of the system. The model is structured 
as a mixed integer prograroming problem vjith integer variables 
taking the values either zero or one. 

The details of models are presented in Chapter 3 

and solution techniques are presented in Chapter 4 which also 
presents the salient features of computer programme. 

Finally in Chapter 5 the model is applied to a large 
scale water resources system, namely, Cauvery river basin 
( South India) and results and discussions are presented in 


Chapter 6. 



Chapter 3 


THE MATHEMATICAL MODELS 


3 . 1 Introduction 

The optimal development of a water resources system is 
formulated in mathematical terms by stating an objective function 
to specify the social goal and constraints to satisfy the * 
different conditions imposed on the model. The locations of 
different projects are assumed to be determined by the topography, 
hydrologi’’ and geological characteristics of the basin. The » 
criterion for the choice among alternative system plans can be 
stated as: select that investment programme and operating policy 
of the system, which -will maximise the present value of net 
benefits, subject to various constraints. The mutual inter- 
action of varionn units of the system at any point in the 
planning horizon is implicitly taken care of in the model. 

The problem under study concerns, v-hich of the projects to 
be built, at vdnat time and hov; to operate the system as a whole. 

The follovring assumptions are made in the analysis; a) In 
investment analysis, a year is taken as time scale of introduc- 
tion of any project into the system (Weiss and Beard, 1970; 
Butcher, Hall and Haim.es, 1969; O'Lagaire and Himmelblau, 1971 ). 
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-o) The demands are taken as deterministic variables. That is, 
the demanci at any point of time is known with certainty. 

c) In power-generation constraints, the head on the turbine is 
taken as constant. 

Three different models are considered for the study. First 
model deals with deterrainistic parameters. Stochastic nature 
of inflows is considered in the sulo— sequent model followed by 
a model v/hich incorTsorates evaporation losses from the reservoir 
surfaces. The models are solved using mathematical programmincf 
techniques presented in Chapter 4. 

3.2 Deterministic Model 

Consider a set of M surface water projects and N groundwater 

projects, each project can be scheduled for construction in any 

one of t time periods (t = 1,2,...T). If denotes the returns 

k-tg 

from project k in period t^ (where k = 1,2,... (M-hlJ)), the probier' 

is to maximize total net benefits from construction of some or 
all of the above projects, subject to various constraints. The 
"renefit from the system is a function of project sizes as well 
as their operating policies, A model representing this objective 
results in a mixed integer prograraming problem, v/here certain 
variables take only zero or one values and other variables take 
real values. 


The problem can be stated mathematically as: 
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Maximize Z, where 


M-fK 


T m+n T M+M T 

^ “ I I A, + I K 1 ~ 1 K 1 gvfCv+. 

o S 


i m-i-n 

I \ I o-,.. 

t =1 -‘V=1 

s 


T M+M 

I \ t, o, 

t„=l ^SK=m+n+l 


(3.2.1) 


SulD j ect to 

1 ) Budgetary constraints 


M+N 

I 

lc=m+n+l ^ s 



< M 


s 


for all t 

s 


(3.2.2) 


2) a) Institutional constraints 
M+N 

I ^v-t- 1 t (3.2.3) 

lc=m+n+l 3 ® 

i.e., number of projects that can be introduced into the syston 
in any year t^will bo at most equal to b, a prespecified imumber. 

T 

b) I g, ^ 1 (3.2.4') 

t =1 k"e k ' s 
s 

v/here k' - a subset of population, is equal to the number of 
projects that can be constructed at a site. This forms mutually 
exclusive project constraint. Only one of k' projects can be 
constructed at a site. For example, at a dam site if three 
heights of dam are possible, then k' talces values 1 to 3 . The 
programme selects only one or none of these k' projects. 
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T 

I Qt--!- < i (for Iv = m+n+1, m+n+2 , ..., M+N) (3.2.5) 

t^=l 

s 

i.e., the project Ic can be built in the time horizon in one of 
the years only, if at all it can be built. 

d) Contingent projects 

(3.2.6) 

S 3 

i.c., the project X is added to the system only if project Y 
is added (Waingartner, 1967') . 

e) Corapound projects 


2^Xt 1 %t 


+ g 


2t, 


(3.2.7) 


Since the value of g is either zero or one, this constraint 
implies that adoption of project X is possible only if the projects 
Y and Z are included (Weingartner, 1967). 


3) Demand Constraints 


a) X. , > Ic = 1,2,..., t = 1,2 (3.2.8) 

ivt j\. L. 

t t 

Q g 

where index t is related to t. by t = remainder of for — ^ O 


t = 2 otherwise. 


is the sum of irrigation, power and municipal and industrial 
and minimum release demands and X_^ is the release. 
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I I ^vt ± 11 ^} 

k t " k t 


(3.2.9) 


For the entire basin total release must be less than the total 
available water in the basin for season. 


4) Phvsical constraints 


a) Minimum storage constraint 


c- ^ -,min 
^Ict - ^kt 


k = 1,2.. .M 


b) Maximurr 


mum storage constrarnl 


(3.2.10) 


< V. Ic = 1,2... M 


(3.2.11) 


c) Conti r 


ntinuicy e^quatron 




R, - for k = 1,2,...M 


kt-1 "let kt 


t =1,2 


(3.2.12) 


d) Canal constraints 


L < X < K for t = 1,2 
m ~ mt — m 


= 1,2... Arcs 


(3.2.13) 


c) I 

' '«* J 


lx.. 


= O for all i 


(3.2.14) 


5) Returns , can be given by the equation 


^ct = “i ^Ict 1. \t + -let 

t=.L 


(3.2,15: 
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6) = 0 or 1 (3.2.16) 

and all other variables are non-negative. 

v/here 

t = tirae index t = 1 - Kharif season 

= 2 - Rabi season 

= time index (t^ = 1,2,...T) years 

index t is related to t bv 

s 

t = reraaindor of t /2 for t /2 O 

S‘ s 

= 2 otherwise 

T = time horizon - years 

k = project index (Ic = 1,2,... m+n+1, ... H-i-N) 

m = nunioer of existing surface reservoir projects, if any 

n = number of existing v/ell-f ield-units, if any 

M = total number of surface reservoirs in the basin 

N = total number of well-field-units in the basin 

M-m = number of nev^ surface reservoirs in the basin 

N-n = number of nevj v/ell-field-units in the basin 

= discount factor 
t 

® 1 

(1+ r ) 

r = discount rate 

A = returns frora project k in period tg 

Ktj 

g^ = decision variable 

J C "C 

= 1 if project k is introduced in period ts 

= O otherwise 
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"kt 


'kt 


X 


kt 


kt 


R. 


kt 


kt 


,rain 

'kt 


X . . 

ij 


' kt 


W. 


kt 


i 

0 ) 

P 

0 ) 

w 


a 


ij 


= capital cost of project k in period tg 

= operation, maintenance and replacement costs of project 
Ic in year tg 

= release from project k in period t 

= storage of reservoir k in period t 

= inflov/ to reservoir I: in period t 

= demand for vrater at project k in period t 

(sum of irrigation demand, power demand, municipal and 

industrial demand and down stream requirements) 

= minimum storage in reservoir k in period t 

= maximum storage capacity of reservoir k 

= flow in the arc a. . 

ej 

= irrigation supply from project k in period t 
= water supply for povxer from project k in period t 
= Municipal and industrial water supply from project k 
in period t 

= return for unit of water (/ Hect-m) supplied to irricatio 

= return for unit of v/atcr supplied for povjer 

= return for unit of vjater supplied for municipal and 

industrial supply 
= arc connecting nodes i and j. 
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3.3 Stochastic Approach 

In the analysis of water resources problem some estimate of 
hydrological, social, technical and economic parameters is 
required. In the previous section - in deterministic model - 
all the inputs and demands are taken to be Icnowii V7ith certainty. 
But these hydrological, economic and other parameters are 
stochastic in nature, Loucks ( 1969) defined uncertain event 
as an event in which even the probabilities of various probable 
outcomes are unknown, and risky events as those events for xdiich 
the probabilities are knovm or can be estimated from the past 
outcomes. Though the demands are also random in nature the 
planners of water resources systems have suJostantial control 
over the outflows. In general outflows can be controlled by 
choosing appropriate types of demand. Therefore, the programming 
technique used in the present analysis incorporates the stochastic 
nature of inflows only. 

3.3.1 Matheraatical programming under uncertainty 

There are various methods of dealing with ■uncertainty. 

These methods are: 

i) sensitivity analysis or parametric prograi'aming 

ii) chance constrained approach 

iii) tvra- stage programming under uncertainty 

iv) expected value approach 

v) minimization of variance approach. 
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The choice of which approach to use depends upon the nimiber 
of factors in the model and computational difficulties encountered. 
The chance constraint method is adopted in the present study for 
the follovring reasons; 

i) By chance constraint method, the stochastic model will 
be reduced to an eoui valent deterministic model, 

ii) Near optimal deterministic models are more easily derived, 

iii) The concept of risk involved is farailiar to water resources 
planners and do not require a sharp break with the tradi- 
tional approaches. 

The formulation of this approach as given by Charnes, Cooper 
and Symonds (1959) can be stated as * 

Max f (C'X) (3.3.1‘) 

Subject to F (AX ^ b) ^ (3,3.2) 

v/here X = decision variable 

C = return coefficients (C' = transpose) 
b = resource value 
P = probability 

= vector of probability. 

In general the decision variables X are selected by a rule 
depending on stochastic variables A, b, and C. Thus X = f(A,b,c) 
and the form of ’f' is often specified in advance - say. linear. 

The linear decision rule (L.D.R.) is used to formulate the 
water resources system managem.ent and design problems as a chance 
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constrained linear prograinming pro?olem. Revelle et al,, (1969) 
used L.D.R. for a single reservoir management model. This 
model was extended to a four reservoir system by Nayalc and 
Arora (1971 ). 

3.3.2 Linear decision rule (L.D.R.) 

The simplest torm of linear decision rule applied to water 
resources system is 

^kt "^kt-l “ \t- (3.3.3) 

where is the release from reservoir k during the period t, 

^kt-1 storage in reservoir k at the end of period t-1 and 

is a decision variable. Equation (3.3.3) states that the 
discharge X, ^ is a function of the storage at the end of the 
previous season since b,,^ is constant. The discharge, hov7ever, 
is not a constant for a given season because the storage at the 
end of the previous season varies. 

The b, , have the follox'/ing characteristics: 

Kt 

i) The release varies inversely v/ith b^,^. In general large 
implies a small release. 

ii) Negative b-, are allowed,- this situation hovrever implies 

j v L. 

that a minimum natural inflow during the season is needed 
to meet the commitment, since the storage 

kt* 


sufficient to insure tlie release X. 
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iii) A negative does not necessarily imply on enpty 

reservoir at the end of the season; if indeed the natural 
inflow to the reservoir were larger than the minimum 
inflov; required to satisfy the release a surplus 

would be generated and stored in the reservoir. 

The system of reservoirs in any basin can be divided into 
two groups, i) independent reservoirs, and ii') dependent reser- 
voirs. In case of independent reservoir generally there will 
not be any reservoirs on the stream upstream of it. The inflox-rs 
to the reservoir is only natural flox^s. In case of dependent 
reservoirs the inflows consist of natural inflov/s plus certain 
portion of releases from upstream reservoirs (fig. 3.l). The 
linear decision rule for each type of reservoir is derived in 
the following sections. 


3.3.3 LDR for independent reservoirs 


In the case of independent reservoir (fig. 3.l) inflox^/s to 
the reservoir are only the natural inflows, the 'net initial' 
storage is equal to their initial storage itself. The continuity 
equation for independent reservoirs is, therefore. 




^Qt-1 + ^Qt 


- X 


Qt 


(3.3.4) 


and LDR is 


y __ Q ^ 

Qt ~ Qt-1 Qt 


(3.3.5) 
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1,2 Independent reservoirs 
3 Dependent reservoirs 


F]G-3-1 SCHEMATIC SKETCH OF INDEPENDENT 
AND DEPENDENT RESERVOIRS 
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= storage in independent reservoir Q in period t 
= inflov/ to reservoir Q in period t 
= release from reservoir Q in period t 
= decision variable of reservoir Q in period t. 

From equations (3.3.4) and (3.3.5) it is obtained that 


S = R + b 
Qt ^Qt --^Q 


( 3 . 3 . S) 


Similarly for time period (t-l). 


S = P 

Qt-1 Qt-1 


'^Qt-1 


(3.3.7) 


i'rom equations (3.3.5) and (3.3.7) the expression for will be 


Qt '‘Qt-l Qt-1 ~ Qt 


(3.3.3) 


3.3.4 LDR for dependent reservoirs 


The ‘net initial storaae' in this reservoir is 


Rt-1 


+ L 


(3.3.9) 


Q u It 


where S^, = storaoe in dependent reservoir at the end of t-l 

R'^ = fraction of the release X_, from upstream reservoirs 

Qt fjt 

Ic that enters the reservoir R. 
up 

LDR is 


^t " ®Rt-l ^ ^ 1 
Q e 1' 


^ ^Qt -^Rt 


(3.3.10) 


Continuity equation 


Rt " Rt-1 


^ ‘'Qt ^Rt ^Rt 


(3.3.11) 
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where z= natural flow into reservoir R in period t 

= release from reservoir R in period t 

From equations (3.3.10) and (3.3.11) the expression for S. 


S - R 4 - 

Rt “ Rt ■ “Rt 


For time period (t-l). 


^Rt-1 ~ "\t-l 


From equation (3.3.10) and (3.3.11) we get. 


-Rt = ^Rt 


I 


"Rt-l ‘ “ '^Rt 

'up 


Rt 

(3.3.12) 

(3.3.13) 

(3.3.14) 


But 


^Qt - 


if' X 
Qt Qt' 


O < lii < 1 

- Qt - 


(3.3.15) 

where = fraction of release that becomes stream flow. 

(3.3.16) 

Osh 
- ' up 

From equation (3.3.8) 


^Rt - ^Rt-1 '^Rt-1 ^ ’^Qt ^Qt ’^Rt 

'^up 


x^^ 

= R_ , "r b 

Qt 

Qt-1 

X 

- R + b 

Rt 

Rt-1 


- b 


Qt 


F (R + b - b ) - b 

Rt-1 ■ Qt Qt-1 Qt-1 Qt Rt 

(3.3.17) 


On simplifying we get, 
X, 


Rt ~ ‘^t-1 ^ Jy ’^Qt ^Qt '"^Rt-l ^^Rt ^ Jy ’^Qt^^Qt-1 "'^Qt^ 

up ■^^''up 

(3.3.18) 


Summarising 

for independent reservoirs 
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Qt 


R 


Qt 


Qt 


Qt = ^'Qt-1 ^^^Qt-1 


for dependent reservoirs 


'Qt 


(3.3.6) 

(3.3.8) 


Rt 


Rt 


-1- b 


Rt 


X 


Rt 


"Rt-l 


(3.3.12) 


'Kt-1 '"Rt Q^T^’^Qt ^‘tOt-l ^ ^^Qt-1 “ '^Qt^ 
Up 


(3.3.18) 


3.4 Stochastic Model 


Objective function 
Maximize Z v/here 


Z 



m+n 


I 

k=l 




I 

lc=m+n+ 


^Ict 

1 


o 




M+N 

I g. 


k=m+n+l ' 




I )^t 

t=l 


s 


m+n 

y 

k=l 



T 

I 


t=l 



M+M 


I 



(3.4.1) 


The budgetary, institutional constraints remain unchanged 
since they contain no stochastic variables. (It is conceivable 
that the outlay in a year could be stochastic; hov/ever, such 
contingencies are not considered here.) 
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These constraints are: 

1 ) Budgetary constraints 


I 

1 v=iTi+ri 4-1 




"s ~s 


fox' all t 

JVU ivL. — ^ 


n 


(3. 


2) Institutional constraints 
h+N 

I £ h for all t^ 

k=m+n+l ' s " 


(3, 


in 1 I = 1 

t=l Jcek^^'^s 
s 


(3, 


T 


11 


Li') I g 


£ 1 for I'C = m+n+1, ... M+N 

t=l " s 


(3. 


iv) ^ for all t. 


(3, 


v) Ig^^ £ ^Zt all tg (3. 

s "" '"s 

3) Maximum storage constraints 

— kt 

4) Canal constraints 

,■ 'l T < n and m (3. 

m — tm — 

ii) Ix.j- Ix-'=0 for all i (3. 

j ij j Ji 

5) Demand constraint 

To provide water for different purposes, the release 
ment should exceed the demand with reliability say Y-, . 

ivt 1 


4.2) 

4.3) 

4. 4) 

4.5) 

4.6) 

4.7) 

4.8) 

4.9) 

4.10) 

comrrd 


In 
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probability notation this constraint is 


P (^rt ^ ^.t) ± Y, 

Similarly ror tne minirnmn storage constraint 

^ ^ \t - 1 Y2 

^//here ^ xs 3.gain irellag--' i i 

^ •— . j. -1- ^ i^Y ^ 


-Y usxng LDR these two constraints (3.4.11) and (3.4.12) 
are converted into eguivalent deterministic constraints. 

3.4.1 Minimum release constraint 

Substituting (3.3.8) in^o (3 4 n) ■>. 

. . U.4.11) rhe release constraint 

or independent reservoir becom'-s 


(3.4.11) 


(3.4,12) 


^ ^V-i + bt-i - t 

1 (Dq, + 

‘V + V - ^Qt-R i (1 - 

V + V - V-1 (1 - t’ 


(3.4.13) 


(3.4.14) 


Substituting (3.3.18) into (3 a n 

. .. 11 .. the release constraint for 


dependent reservoir becomes 


Rt-1 ^ ■‘^Rt-1 "Rt ^ (Rq^ . + b T - b ) 

Oe]c Qt-1 Qt-1 "^ot^ 


■"Rt ^ - a 


(3. 4.15) 
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i.a., P{R 


■Rt-1 — 


t ^Rt ^Rt-1, ~ I ^r,4- -, -h )\ 

up 

-^'l 


^ "’■*^R(DRf*Rt+]-Rt-l ~ ^ ^ (R '-h ’ A ^ 


Fri ■f + n ■''' 

R Rt -"Rt -~^F 


Oek '-Qt-l-^Qt_l--Opt^^ > Y, 

up 


i!) Ct> 
VK, 


Qek Qt -^Qt-i 


Qt-1 “ 


1(1 - Y •) 


- C. x\*( 


Rt-1 


I 


Qt %t-i + -V_1 - F Ri-Yj3 


(3. 4.16) 


3.4.2 M 


inim-uni storage constr 


raint 


Substituting (3.3.6) into (3.4.12) the 


constraint for inder 


minimum storage 


pendent reservoi: 


P (R 


, > e^nin,, ^ 

Qt ^ -^ot - ^O-r ^ 1 Y. 


(3. 4.17) 


P «Qt i - >=.30 ) 1 V. 

" - k - >=ht> 1 V, 


oiTiin , _1 

SqO - >=Qt i P 


(1 - vg 


(3.4.18) 


Similarly for dependent reservoir the minimum storage constraint, 
from (3.3.12) and (3.4.12) will be 


P (R^, + b, > > Y 

j3t Rt — Rt — '2 


(3. 4.19) 


=Rt - ^ (1 - Ig) 


(3.4.20) 
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where (^1 is a solution of 




P(r) = F(R_<r-) _ Y 

(3.4,21) 

(r) is a .nown ^aantity that can he obtained fnc the cupula- 
rive distribution function of inflows. 

Thus the chance constraints are replaned 
n , . . ir' ea D_i t.ne equivalent 

deterministic constraints. 


1.5 Chance Constraint Formulation with Evaporation losses 

in and and serai-arid regions, evaporation losses play 

important role in planning and designing of water resources 

systems. Generally the evaporation losses an o 

losses are expressed as 

percentage of storage in any time period in thr 

in the present analy= 

the expression of the type S, = 5 (3 ^3 . 

■■^t ht '■ bt ^ ]rt-i^ used 
(Chow Windsor, 1972), ^here K j c, 

X. ere is the evaporation losses in 

period t ror reservoir k, and 6 ^,^^ evaporation coefficient 
As Harbech (1956) suggested the depth of evaporation loss in any 
IS -clouiiely constant from year to year, i.o., the valu 
of tt -mains constant from year to year. Cho evaporation 
losses can be expressed as a linear function of surface area. 

But at any specific reservoir site there is a definite relation- 
ship between surface area and storage volume. Hence the evapora. 
tion loss expression is expressed as storage volume and approxl- 
mated by a straight line with positive slope. 

In the present analysis, a projected storage based on the 
actual storage at the beginning of the period and on the anticios 
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evaporation loss during the period is taken for formulation 


the model. 


expression for evaporation losses is 




(3.5.1) 


The linear decision rule can be extressed 


kt ict-1 


(3.S2) 


v/nere 


kt-1 


13 projected storage taking evaporation into 


consideration. For the linear decision rule to be an effective 
computational device, is expressed as a function of the 

storage (Rovello & Kirby, 1970) . 


S, . = , S. , T 

ku-l ku Kt-1 


(3.5.3) 


where a 


1 - 6 . 


(3.5.4) 


The continuity equation for the reservoir is 


Sv., = S, , , p R, ^ _ X, 

kt kt-1 let K 


(3.5.5) 


Substituting (3.5.1) and (3.5.2) into (3.5.5) we get 


R 4 * 

kt ' kt 


(3.5.6) 


= {R, . - b, , ) Gt 

kt let ku 


(3.5.7) 


where , = 1 + , 

kt kt 


(3.5.8) 


Similarly, + b,,^_^) 


(3.5.9) 


and f rom (3.5.3) and (3.5.9) 


kt-1 k 


“ ’^kt ^ ^lct-1 ^kt-1^ 


(3.5.10) 
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kt 




{R-, + b, } - b^^. 

ivr-1 Kt-l i<^t: 


(3.5.11) 


3.6 Model v/ith Evaporation Losses 

The objective function (3.4.1) and constraints (3.4.2)- 
(3.4.10) remain the same. The minimum release constraint and 
minimura storage constraints become: 

Minimum release constraint 




SUiOstituting for X,,j_ from (3.5.11) 

it L. 


P { 


kt 


( D 

'• * ^1,. j 


®kt-l 


+ b 


j _ b > ^ Y^ 

ct-1^ kt'' Icc 1 


e 


P { R. 


kt-1 


kt-i i'hj 


kt 


kt-1 

^kt 


> > \ 


e. 


p { R, 


kt-1 


i (Eivt + 

ic-c riT . 


'kt-1 ^ - "^1 


'r ^°kt ^ '"^kt^ 




lct-1 
^1. 




kt 


]ct-l - 


< (1 - yJ 


e. 


(D + b ) ■■ - - - - - b < F~^ (1 - "'^1 ^ 
^ kt ^ kt^ nv^ kt-1 - 1 


(3.6.1) 


(3.6.2) 


(3.6.3) 


Minimum storage constraint 


P (S^ 


kt 


_min' 

'kt 


> s;“r^) 


> Y, 


(3.6.4) 


Substituting for from (3.5.6),& (3.6.4) the follovring 


expressions are obtained; 
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The suJoscripts in these expressions are suitably modified 
for independent and dependent reservoirs as explained earlier in 
section 3.4. That is, for independent reservoirs the subscript 
k should be replaced by Q, and for dependent reservoirs it should 
be replaced by R. 



Chapter 4 


METHODOLOGY 

BRANCH AND BOUl-ID AND OUT-OF-KILTER ALGORITHMS 


4.1 Introd-uction 

The methods generally employed in water resources systems 
analysis are mathem.atical programming m.ethods and simulation 
techniques. Simulation is a process that duplicates the 
essence of a system witiiout attaining reality itself (Tocher, 
1963'). A system is modelled either on analog or digital or on 
both. The response of the system to any input data can be 
predicted. However, the simulation model is not self-optim.izing. 
Simulation is an excellent tool in modelNnc complex systems 
but has less utility where the number of variables are large 
as in case of v/ater resources systems. Mathematical programming 
methods are self-optiraizing. There are many programming methods 
which are used in the water resources systerrs analysis . 

In Chapter 3, the water resources syston model is 
developed as a combination of investment policy problem and 
operating policy problem v-hich are solved iteratively. 
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The problem is of zero-one (0-1 ) mixed integer programming 
type where the system elements are either selected or rejected 
and operation policy giving the releases to meet the various 
demands is obtained. Consequently, all the feasible solutions 
contain a mixture of integer and continuous variables, and 
integer variables are restricted to the values zero or one. 

In solving the model it is required to find (i) which of the 
projects should be built, at what time and to what scale, and 
(ii) the amount of releases such that the objective function 
is maximized. 

Mixed integer programming method requires that the objective 
function be separable with respect to integer and continuous 
variables, i.e., the two sets of variables must be capable of 
being ' linearly sumraed (Hadley, 1964). But in the model presented 
in Chapter 3 there is an interaction between the variables 

in the objective function and also in some of the constraints. 
Therefore, a nev; method of approach is sought for solving the 
model . 

The problem is decomposed into two units as a combinatorial 
problem (a) investment problem, and (b) operation policy problem. 
The investment problem is solved by branch and bound technique. 
The operation policy problem is solved by net-worlc analysis 
using Fullcerson's out-of-Kilter algorithm. The connection betv/oe 
these two problems is due to the fact that the investment policy 
problem is solved by obtaining optimal solution of operating 
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policy problem and the operating policy problem is dependent 
on the number of projects obtained by the investment policy 
problera {Figure 4.1; . 


i -i-nvestaent policy problem ! 

r * 

I Maid-mize net benefits of 

! all feasible proiects sub- 

I 

: joct to budgetary;- & insti- j 
• tutional constraints I 


1 

; Kumber of 
'projects in 


i the system 

r 


Optimum 
operating 
policy of 
the system 



Operating policy problem 

Maximize the returns from 
operating the selected 
system from investment policy 


! problCTi subject to demand 



■' suppl y a nd physical constrain ts 



Fig. 4.1 - Connection betueen Investment policy and 
Operating policy problems. 

4.2 Investment Policy Problem 

This problem falls into a category of integer programming 
methods with -variables taking only tv/o values, either zero or one. 
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The problem is solved by branch & bound technique (Little et al., 
1963). Branch and bound is a tree searching technique. The 
common features of this technique are: 

(1) They are easy to understand, 

(2) They are better suited for computer application than 
cutting plane techniques. 

Branch and bound technique is structured as a search 
technique in the space of feasible solutions. The feasible 
solution space is repeatedly partitioned into smaller and 
smaller sub-sets, and the bound is calculated for each sub-set. 
After each partitioning, those sub-sets v/ith a bound less than 
that of the known feasible solution are excluded from further 
partitioning. The partitioning continues till the feasible 
solution is found such that its return is crreater than the 
returns of any sub-set. 

For zero-one integer prograraming with p variables there are 
2^ solutions. These 2^ solutions are partitioned into n+1 sub- 
sets. The s suls-set (s = 0,1, 2.... n) contains all solutions 
with s variables being 1 and (p-s) variables being zero. The 
zeroth sub-set contains one solution yj=0 (J = l,2,3,...p). 

The first sub-set contains (^) solutions: Y^=l and 7^=0, (j^^l); 

y„=l, y.=0, and so on. In general, s^^ sub-set has (^) 

z J S 

solutions. In testing the feasibility of solutions, the values 
of variables y^. are substituted into the constraint set directly. 
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The process starts with the highest node; and then follows 
an intermediate successor to the highest node. A node is said 
to be a terminating node if any of its successors nodes are not 
tested. 

The investraant policy problem is represented as in Figure 

4.2. 

The operation policy problem is solved for the original 
systen as v;ell as the system with addition of a project in any 
year t. If, and only if, the solution of latter operation 
policy problem gives a higher return than the operation policy 
of former problem, the new project is selected in the year t. 

The procedure is repeated for each subsequent year. The termina- 
tion of branching occurs at T (tirae-horizon) or v/hen all the 
projects have been introduced into the system. The optimal 
solution is found when all the feasilcle nodes have been searched. 

4.3 Operation Policy Problem 

From the investiaent polic^p if it is decided to introduce 
a nev? project into the system, at the beginning of year t, its 
feasibility will be tested x#ith respect to its operational policy. 
The operation policy problem is solved as a net-work problem. 
Fulkersons out-of-Kilter algorithm is used in the analysis. 

The follov/ing gives a brief description of the method. 

A net-work consists of a set of nodes and a set of arcs 
connecting the nodes. Figure 4.3. 
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suoispap ajqtssod 


Time in years 
f m + n Existing 
F Total number of 
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Fig, 4.3 - A typical net-work 

If the arc has a direction then it is called a directed arc 
In water resources systems the directed net-worlc analysis is 
used. There are tv/o special nodes in the net-worlc, namely, S_ 
(source) and S-,. (sink). Let S. and S. be the nodes and a. . be 
the arc joining S. to S.. The sequence of nodes and arcs S„, 

1 J K. 

^Rl' ^1' ^12' ^2 *•*' is called a chain leading from to 

S-,,. If S„ coincides with S, , then it is called a directed cycle 

Jv K iv 

Associated with each arc there is an upper flox"/ capacity H. . 

Ij 

and a lov/ flow capacity L._. and also the return b. . from passing 

— J 
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one unit of water through the arc. Let the flow in the arc a. . 

ij 

be X. .. Then the problori is to find the values of x. . such that 

ij 

the returns will be maximijm v/ithout violating the constraints. 
Maximize Z 

v/herc Z = ^b.^. x . (4.3.1') 


Subject to X. . < K. . 

^ Ij - XJ 

for 

each a . . 
ij 

(4.3.2) 

-X , . < -L . . 

- IJ 

for 

each a . . 

( 4-. 3 . 3 ) 


11 

I ^ 

X . . 

= 0 for each i 

(4.3.4) 

j j 




X . . ^ C for each a. . (4.3.4) 

From duality theory, there is a dual variable associated with 
each primal constraint and a dual constraint associated v;ith each- 

primal variable. Let o_. denote the dual variable associated 

tn. 

with the i " primal conservation of flov/ equation (4.3.4). Let 

y^j denote the dual varialole corresponding to upper bound 

constraint (4.3.2) and Z. . be the dual variable associated with 

ij 

lov/er bound constraint (4.3.3). 

The dual problem can be stated as optimize where 



V 

L 

j 


H, 


'iJ 


I I 

i j 


L. . Z. . over all arcs 
ij ij 

(4.3.6) 


+ y. . - z. . > b. . for all i,‘j‘' ’• ' (4*.§P/5 

1 J ij - ij 


Subject to 
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y . . > 0 

2 > Q for all i £i j (4.3.8) 

ij - 


and the variables are xxnrestricted in sign since these dual 
variables are associated with quality constraints in the primal 
problem. 

At the optimal value, the value of primal and dual objec- 
tive functions are equal ( Z = Z' ) . The relationship between 
the primal and dual variables is known as complimentary slacjcness 
conditions (Hadley, 1964') . These conditions are 

if y . . > O, it follows that x..=H., (4.3.9) 

ij iJ 

if z. . > O, it follows that x. . = L . (4.3.10) 

-L X J_ 

and if a - a + y - z . = b it follows that b . £x. .£H. . 

1 J Ij IJ Ij IJIJIJ 

(4.3.11) 

Fulkerson (1961) developed the out-of-Kilter algorithm (O.K.A. ) 
which takes the advantage of above relationship and examines a 
relatively small sub-set of the dual variables in a search to 
get optimal value of the objective function. The efficiency is 
achieved by defining the quantities. 


C. . = a. - a. - b. . (4.3.12) 

ij 1 J 11 

y. . = Max (0, - C. .) (4.3.13) 

^ij IJ 

Z. . = Max (0, C. .) 
ij iJ 


(4.3.14) 
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such that y. . and z. . satisfy (4.3.7) and (4.3.8). With these 
definitions, the complementary slackness conditions may be 
written in the follovring form. 


if 

c. . < 

IJ 

O, 

X . . = M . . 
ij 13 

(4.3.15) 

if 

c. . > 

IJ 


X. . = L. . 

13 13 

(4.3.16) 

and if 

c. . = 

0, 

L . . 

13 ''i3 ■ ij 

(4.3.17) 


These three conditions along with continuity equation 
condition (4.3.4) determine the sufficient conditions of opti- 
mality. 

Those arcs v;hich do not satisfy at least one of the above 
conditions are known as out-of -Kilter arcs. The general method 
of solving net-work problem is given by Fulkerson (1962) . The 
O.K.A. seelrs systematically to direct the flov/s and assign the 
values to the o- variables so that each arc satisfies the opti- 
mality conditions, i.e., each arc is in in-Kilter condition. 

The variable o- ^ may be considered as the price associated 
v/ith a node i and represents the total price of unit flov; 

from node i to node j . 

With the above optimality conditions, different states of 
arc stated as in the table 4.1. 
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labelling the nodes to indicate the direction that the flow in 
an arc may be increased. Once the flow augmenting path is 
found, tne flow is increased to its maximum capacity along the 
path and all the labels on the nodes are erased. Then assign 
new laJoels to nodes based on nev/ flovjs. Thus this algorithra 
consists of tv/o steps. Step 1 is the labelling procedure and 
step 2 is the implementation of flow change 

■Step li Labellinc- Proced'ure 

Every node is always in one of the following three states: 
a') Labelled and scanned 
b') Labelle) and unscanned 
c) Uni abell ed . 

A label for a node always has tv7o parts, first part is 

J 

the index of node from v/hich the flow can be sent to S,.; 
and second part is a number v/hich indicates the maximum amount 
of flow’’ that can be sent to Sj v/ithout violating capacity 
constraints . 

A node is labelledt and scanned if it has label and all its 
neighbouring nodes are inspected. 

A node is labelled and unscanned if it has a label and not 
all its neighbouring nodes have been scanned. 

A node is unlabeHed if it has no label. 

Initially all nodes are unlabelled. Select a node and 
assign the label of the form (i'*', e (j)) or (i , e (j)) where 
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states that the flow can be sent from i to j and '^{j) 
indicates the amount of flov/ that can be sent. 

To all nciahbouring nodes 3 of 3., which are unlabeled, 

m j 

and for v/hich O < x. < H. , assign the label (j'^, e (m) ) where 

— im — jm 

^ (ra) is equal to min {e(j), H -x . } and to all neighbouring 

jm jm 

nodes S vrhich are unlabelled and for vjhich x . >• O, assign the 

rn [H J 

label { j , e(ra) } where e (m) = min { e(j), x . The positive 

m J 

( 4 -ve) and negative (-ve) signs in the label indicate, how the 

flow should be changed in the step two. Novj v/ith all neighbouring 

nodes of S, which are scanned, 3. is considered to be labelled 
J J 

and scanned. All nodes 3 are nov? labellcsd and unscanned. Continue 

m 

to assign labels to neighbours of labelled and unscanned nodes, 
until either S^. (sinlc node) is labelled or no more labels can be 
assigned, and 3.,, is unlabelled. If 3^, cannot be labelk?d., no flow- 
augmenting paths exist and hence the flow is maximum. If 3, is 
labelled, a flow augmenting path can be found using step two. 


Step 2: 


Assume that S-,, is labelled (m”^, e (k) ) . Replace x^^^, by 

X , -i- ^ (k) and turn to 3 . If S is labelled (j”*", e (m) ) , replace 

mk m m 

X. by X. -i- £(m), and turn to 3 . . If S^, is labelled (j~, e (m) ) , 
jm - jm j 

replace x . by x . - e (m), and turn to 3.. Continue this procos: 
- mj ^ mj J 

until is reached. Erase the labels on all nodes and go back 
to step one- 
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Steps one and two are iterated until the maximum flov/ condi- 
tion is obtained. Out-of-Kilter algorithm steps are given in 
the flow ciagrarTi/ Figure 4,4. 

A computer prcgrarome is available (Himmelblau, 197 3). It 
is basically in the original form as v/ritten by Texas-Water- 
Development- Board and apparently has not undergone testing on 
large scale problems. Several errors in the prograiiime have been 
detected by the writer during this v/orlc and several subroutines 
are introduced to suit the present large scale river basin model. 

The model is solved using IBM 7044/1401 system. 



Initial flow 
Assiqments 


Cjj =a j-aj -b ,j 


rase a 


W h 6 1 h 

all arcs are in 

Kilter cheek the 
conditions^^ 


Change the flow 
along the path 


Labeling & 

flow augmentation 
process 

— 

Does^\^ 


The flow augmentation 
path is obtain^ 


Change node ; 
prices and : 
label the nodes 


Whether^ 
ll arcsiri Kilte 


Optimal solutio 


I 


Fie- 4'4 FLOW DIAGRAM FOR OUT-OF-KILT 









Chapter 5 


MODEL APPLICATIOM TO ^ 
CAUVERY RIVER BASIN - (KARNATAKA '-INDIA) 


5 . 1 Introduction 

The different models presented in the Chapter 3 are applied 
to the development of water resources in the Cauvery basin. 
Economy of the region depends almost entirely upon agriculture 
and most of the v/ater is reguired for irrigation piorposes. 

The physical and hydrological characteristics of study area 
are briefly presented first. Then the mathematical model repre- 
senting hydrological, physical, economic and institutional para- 
meters is presented and solutions are obtained. 

5 . 2 Physical Description of the Model 

The study area is that part of Cauvery basin vfhich lies in 
the Karnatalca state. It comprises of upstream and middle reaches 
of the Cauvery river (fig. 5.l) 

The drainage area of the basin under study is approximately 
3 6,000 sq.]<m.. Entire Cauvery basin spreads betv/een the longitude 
of 75° 30 'E & 79° 45°E and latitudes of lO^S *K & 13°30'K. 


*Name of the state was changed from Mysore to Karnataka, 1.11.197 
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Located in peninsular India, the basin covers the areas 
in the states of Karnataka, Kerala and Tamilnadu. 

The state-'vise distribution of drainage area is given in 
the table 5.1 (CFPC, 1972) 


Table 5 . 1 

Area of the State and drainage basin 


State 

Area of the 
state (sq.km. ) 

Area of drain- 
age basin 
(sOa lora. ) 

Remar }cs 

Karnataka 

1,91,733 

36, 240 

Drainage basin 
includes vjcstern 
ghats area v/herv'. 
the river rises 
and the plateau 
of Karnatalca 

Kerala 

38, 864 

2,930 

Drainage area 
consists of ghat 
areas 

Tamil Nadu 

1,30, 070 

48, 730 

Drainage basin 
includes middle 
reaches and 
del ta area 


Total area of drainage basin 87,900 sq. 


The area that lies in Karnatak.a state is considered for 
analysis. This study area is divided physiographically into 
western ghats region and plateau of Karnataka. The v/estern 
ghat area is mountainous and covered with thick vegetation. 
River Cauvery takes birth at Talakaveri in Coorg district of 
ghat areas (Figu 5.l) . The location vThere it rises is roughly 
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25°24'N longitude and 74°34'E latitude at an elevation of 1341 m 
above M.S-L. in Brarahagiri range of hills. The river has steep 
banks and many ideal sites are available for impounding water. 

The plateau of Karnataica ’/jith average elevation of 750 m, 
slopes gently tov/ards S-llE direction. 

There are a niamber of tributaries which join the river along 
its course from both sides of the river (Fig, 5.1'). The 
important tributaries from the left side are Harangi, Hemavathy, 
Shimsha and Arkavathi. The tributaries which join the river 
from right side (south side) are Lalcshmana-Theertha, Kabini, 
Suvarnavathy and Uduthore-halla. 

The catchment area of Harangi river lies in ghat area 
v/hich receives heavi’' rain fall by south-v/est monsoon (June to 
October). It drains an area of 540 sq. len. and joins the river 
Cauvery at Kudige in Coorg district. 

Hemavathy river, with its tributary Yagachi, is a major 
tributaey v/hich joins the river from the left side. The total 
drainage area of this tributary is 5200 sq.km. 

The other two tributaries which join the river from left 
are Shimsha and Arkavathi. These tvjo rivers flov7 in the plateau 
of Karnatalca- which is a rain-shadov/ area of the x-zestern ghats. 

Lakshmana-Theertha river joins the river Cauvery from south 
and drains the areas in Coorg and Ilysore districts. 



63 


River Kabini along with its numerous tributaries is a major 
river which joins Cauvery river from the south. It rises in the 
ghats of Kerala state and drains a total area of 6690 sq.km, 
both in Kerala and Karnatal'a states. It traverses a length of 
210 Jem before joining the Cauvery river at Tirumakudala Karasipura, 
in Mysore district. There are a number of sites available on this 
river for conservation and control of water. 

Suvarnevathy and Uduthor e-hall a are other tv/o tributaries 
which join Cauvery from right. These two tributaries flovj- in 
Karnatalca plateau region and cariry moderate runoff. 

The upper reaches of the basin are steep and river flov?s n 
narrow gorges. At Sreerangapatnam (d/s of the city of Mysore) 
the river divides into two branches and forms an island where 
the famous bird sanctuary Rangana-Tittu is situated. The river 
flox*7s in eastwardly direction in the plains of Karnataka. At 
Shivasamudram a fevr Icilometers upstream of the confluence of 
Shimsha^, the river divides into two branches and falls through a 
height of more than 91 m in a series of falls and rapids. The 
two major falls are Gagana Chu]cki and Bhara Chukki. ’ There is an 
hydroelectric station at this fall; and it v/as constructed in 
1902. The two branches of river join dov/n-stream of the falls 
and flov 7 in a narrov7 gorge. At one reach the width of the 
river is so narrov7, it is naraed Mekedatu which literally means that 
a goat can leap across. After this point the river flows in- 
eastward direction and receives Arkavathi from left and Uduthore- 


halla from the' right. 
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xRiver Cauvery runs a total length of 320 Ion in Karnataica 
state and for auout 64 jnri it forms the boundary between Karnatalia 

and Tamil Nadu States. 

5.3 Climate 

The climate of the basin area under study is essentially 
uropical monsoon, type. The year can be divided into tv/o seasons 
as per agricultural practices - Kharif and Rabi seasons. Kharif 
season spans from month of May to month of October. It includes 
not months foliovred by south-west monsoon period. Ghat areas 
and adjoining upper reaches of the basin get major portion of 
the rainfall in this season. 

Plateau of Karnataka receives moderate rain fall during this 
season. 

Rabi season comprises of the months of November through 
April. During this season north-east monsoon gives moderate to 
heavy rain fall in the plateau of Karnataka in the months of 
October and November. 

The nonual annual rainfall distribution is shov/n in Table 5.2. 

Table 5.2 

Rain fall distrileution 
Month /o of rain fall Remarks 

_ ^ zyy 3 


1 

January 

0.30 

2 

February 

0.40 

3 

March 

0.90 

4 

April 

4.50 


Rabi Season 



65 


Ta-ole 5.2 (...Contd.) 



1 

2 

3 

5 

May 

3. 30 


6 

June 

14.60 


7 

July 

2 5 . 5 C 


8 

August 

16.10 

I'htarif season 

9 

September 

11.90 


10 

October 

1 0. 10 


11 

November 

5.20 

Rabi season 

12 

Dec amber 

1.00 



Evaporation s This is an important parameter in water resources 
planning and management studies. Unfortunately ver^.'' little 
information is available about evaporation losses. Hov/ever, the 
evaporation losses are estimated using meteorological data. The 
average daily evaporation losses are given in Table 5.3. 


Table 5.3 

Average evaporation losses 


6 



Month 

Average daily 
evaporation loss 
in mm 

Percentage 
of loss 


1 

2 

3 

1 

J anuary 

5.50 

7.10 

2 

February 

7.00 

3.90 

3 

March 

7.60 

9.60 

4 

April 

7.10 

9.00 



66 


Table 5.3 (...Contd.) 



1 

2 

3 

5 

May 

7.00 

3.90 

6 

June 

8.30 

10.50 

7 

July 

7.10 

9.00 

O 

August 

3.80 

8.60 

o 

September 

7.30 

9. 20 

10 

October 

5.80 

7.40 

11 

November 

4.70 

6. 00 

12 

December 

4.60 

5.80 


5.4 Geology and Soils of the Basin 

Gneisses cover the major portion of the basin. Granite 
occurs as intrusions in the gneisses shov/ing characteristic 
structure of plutonic rocks. The Arkavathi sub-basin contains 
granite formation which runs from Tumlcur to Sivasamudram. 

The study area is covered with red soil- The loamy structure 
of red soil make it suitable for growing a large variety of crops. 


5.5 Crops 

Rice is the widely grovm irrigated crop follovred by sugarcane. 
Other important crops are ragi, groundnut and some garden crops. 
Who.le of arable area can be brought under irrigation v/ith ass\:r:.i 
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supply of water. The entire basin has excellent drainage facility. 
Therefore, irrigation drainage presents no problon. 


5.6 Water Resources 

The most comraonly used geographic units for water resources 
planning and development is the river basin or a closely related 
group of basins, v/hich drain to a common point. With such a 
nydrologic complex the v/ater supplies are connected. Within the 
rnajor area under consideration many streams and stream systems 
malce up the smaller hydrologic units which lead themselves to 
analysis as individual units. The entire Cauvery basin is 
divided into sixteen hydrological units. Out of these sixteen, 
first seven lie in the study area (Karnataka part of the basin) . 
Table 5.4 and Figure 5,2 give the details of these study units. 


Study unit 


C 


1 



C 


6 



Table 5.4 

Details of Hydr^ical Units 
Area covered 

Upper reaches of Cauvery, Harangi river, 
Hemavathy river system and Laksbmana Theertha 
river 


Kabini river with all its tributaries 
Suvarnavathy and Gundal rivers 
Shimsha river system 
Arkavathi river system 

Area along the course of the river which ir 
not covered by other units 

Uduthore-hal la basin 
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Discharge sites in Cauvery basin 
(Fig. 5.2) 


31 • 
No . 

Gauge station 

Karae of the Hydraulic 
river or study 

tributar^^ unit 

Year from 
vjhich 
gauging 
is done 

Method 

of 

gauginc 

1 

2 

3 4 

5 

6 

1 

Kushal nagar 

Cauvery 

1967 

Qi 

2 

Chunchanalcatte 

Cauvery’’ C^ 

1916 

CM 


Kr i shn ar a j a s ag ar 

Cauvery 

1934 


4 

lannur 

Cauvery 

1971 

Q! 

5 

Dhanagsre 

Cauvery’’ 

1922 


■Q 

T. Warasipur 

Cauvery C^ 

1966 

’CM 

7 

Sathyagala bridg 

e site Cauvery 

1970 


o 

Shivasamudraiu Anicut Cauverg.'' 

1939 

MI 

15 

Karnbibane 

Chicklihole 

19S6 

WA 

16 

Hudgur 

Harangi 

1964 

Wi 

17 

Kudige bridge si 

te Harangi 

1971 

m 

13 

Sakaleshpura 

Hemavathi 

1971 

ai 

19 

Ch i]c] caby ad ag er e 

Yagachi 

1965 

SA 

20 

Belur bridge 

Yagachi 

1968 

04 

21 

Madaghatta 

Votehcle C^ 

1968 

. 

22 

Chennanahalli 

Votehole 

1965 

SA 

23 

Vidyapeeta 

’Tagachi 

1971 

CM 

2 4 

Gorur 

Hemavathi 

1966 

SA 

25 

Srirana Devaru 1 - 

‘.nicut Heraavathi 

1950 

T''T7 ^ 

26 

AkJcihebbal 

Hemavathi 

1916 

C*/' 

27 

Karihar 

Lalcshraanathirtha 

1965 

S':- ' 

28 

Unduwadi 

Lalcshmanathirtha 

1916 

C : 



Table 5. 4 
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(. . .Contd. ) 


29 

orinivasa Agreaiar'a 

Pur a 

H.D. Kote 

Wugu dam site 

Lo kapavani 

Kabini 

Taralca 

Nugu 

4 

5 

6 

33 

3 4 

35 


1967 

1956 

1965 

1961 

SA 

a-i 

Q-I 

RM 

36 

Kul 1 3.11.0.1 9-11101111. 

Kabini 

^2 

1939 

01 

37 

Hullahalli bridge 

KaJoini 


1970 

ai 

38 

Nan j angucl 

Kabini 


1966 

S-A. 

39 

T* I’cirasipur 

Kabini 


1966 

a: 

41 

Kbllegal 

Suvarnavathy 

^3 

1965 

CM 

42 

Marconahally 

Shirasha 


1940 

Rl-i 

43 

Muddanahalli 

Shims ha 


1967 

SA 

44 

Icanva Reservoir 

ICanva 

b 

1946 

Rr? 

45 

T. K. Halli weir site 

Shimsha 


1942 

T,T-,,T 

46 

T,IC. Halli bridge site 

Shimsha 


1970 

G1 

47 

Charaaraj as agar 

Arkavathy 


1939 

RK 

43 

Manchanabele site 

Aricavathy 


1965 

SA 

49 

Kanakapura bridge site 

Arlcavathy 

“5 

1970 

Q-1 

50 

Arobele 

Arkavathy 


1965 

SA 
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Table 5.4- ( . , . Contd. ) 


1 

2 

3 

4 

5 

6 

51 

Chengawadi 

Thattehalla 


1961 

ai 

52 

Rarapur bridge 

site Uduithorehalla 

s 

1970 

Q1 

53 

Manuganahall i 

Doddihalla 


1961 

SA 

5 4 

Uggyam 

Hinnattuhalla 


1966 

SA 


Observations given under column 6 indicate the method of gauging 
as; 

O'! Current meter 
WM Weir calculation for flow 
RM Reservoir operation 
3A Slope area method 
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The locations of rain gauges and stream guage stations are 
shown in iig, 5.2. The rain fall and run off records are available 
for most of the stations. 

5.6.1 Ground water resources 

ihe geological rorrnulation is not favouraJole for ground 
vjater storage. However, some of the valley portions are filled 
V7ith wash derived frora the adjacent highlands- The area along 
the banlvs of some streams form a good aquifer for storage. 

No systematic study has been done for the assessment of 
ground v/ater potential in the basin. Raghava Rao et al., (1969) 
located some poc''.ets of ground v/ater areas. These are in tlie 
hydrologic study units, C^, in Kemavathy sub-basin, and 
in Shimsha and Ar]'avathi basins; and in Uduthore-halla basins. 

The groud v>rater potential of these regions is given in table 5.5. 

Table 5.5 

Ground Water potential 

Area Ground water potential in 100 Hect-rn 

- Hemavathy basin 90.00 units 

Shimsha and 226.00 units 

ArJcavathi basins 


- Uduthorc-halla 


136.00 units 
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Generally open dug wells are used for exploitation of ground 
water. Several of these v/ells put together is considered as a 
unit. a.he locations of the units are shown in flovr diagram 
(fig. 5 . 3 ) . 


5 . 7 Mathematical Model 


i'lov: diagram (fig. 5.3) shoves the main river with its 
tributaries/ the location of various projects, irrigation 
districts, hydiro elec trie power stations and municipal and 
industrial deraand locations. There are 100 nodes indicating 
junctions and diversions. The model is studied for tv 7 o seasons, 
both Kharif and Rabi. The node numbers of figure 5.3 refer to 
those of iCharif season. The second season's (Ralai season) nodes 
are obtained by adding an arbitrary number say 150 to the nodes 
of the first season (Idiarif ) . There are 437 arcs in the complete 
model (for both seasons) . The flow diagram is summarized in 

* 

table 5.5 to indicate the flox-Js to and from the nodes. The 
number in a row with a positive (+vc) sign indicates the direc- 
tion of flov/ to the node and the number v/ith negative (-ve) sign 


indicates the direction of flov/ out of the node. Tor instance, 
at node 3 (column l) the direction of flow is from the nodes 1 
and 2 vjhich have positive signs and outflovj direction is to the 
node 7, which has negative sign. There is no power or irrigation- 
diversions from this node 3. 


There arc two super nodes namely source node (299') and sinh 
node ( 300 )-. All inflov;s to the system and original reservoir levels 
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Table 5.6 

Flow Diagram Details 
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TaJDle 5.6 ( . . . Contd. ) 
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originate from the source node. All exit flows and final 
reservoir levels are removed to the sink node. An artificial 
is introduced connecting sinlc node to source node to 
complete the circulation network. Figure 5.3 gives the details 
of flov7 inputs and outputs of the system. 

5.3 Economic Characteristics of tiie Systen 

Costs and benefits of the system are required to formulate 
the mathematical programming model. The capital cost, and the 
operations, maintenance and replacsnent costs are obtained from 
the project reports. The costs of vjell projects are estimiated 
after getting relevant data from individual well owners and 
officials of Central Ground VJater Board. 

The benefits from irrigation are estimated using land 
revenue values since v/ater rates are collected in the form of 
land revenues. The returns from pov/er and municipal and indus- 
trial water supplies are estim.ated using the data provided by 
Karnataka Electricity Board and Bangalore Water Supply & Sewarage 
Boards (Appendix A). Nayak and Arora (1971 ) reported that it is 
enough to Icnow the relative economic values rather than actual 
values of the projects for scheduling and sequencing problems. 

The major uses of water in Cauvery river basin are irrigatic 
hydroelectric power and municipal and industrial water supplies. 
Other uses like recreation and ecology are implicitly taken 
into consideration by providing minimum pool level in some 
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reservoirs. -irrigation is a major sector which requires large 
percentage of v/ater availatle in the basin (nearly 70% or more) . 

5.9 Existing Facilities 

There are 12 existing reservoirs which supply v^ater for 
irrigation/ municipal and industrial supplies and for pov/er 
generation. Table 5.7 gives the details of these projects 
(figs. 5.3 and 5.4). 


Table 5.7 
Existing projects 


SI. No. 

Node Mo. 

Pro j ect 

Capacity 
(100 Hect-ra) 

Purpose 

1 

17 

Hemavathy 

960.00 

Irrigation 

2 

21 

Krishnaraj asagar 

1140. 00 

Irrigation, pon 
recreation 

3 

31 

Kebballa 

11.00 

Irrigation 

4 

35 

Kabini 

541.00 

Irrigation & 
pov/er 

5 

42 

Nugu 

122.00 

Irrigation 

6 

58 

Chil'hahole 

10.00 

Irrigation 

7 

63 

Marconahal 1 i 

64.00 

Irrigation 

8 

66 

Mangala 

7.00 

Irrigation 

9 

71 

ICanva 

18-00 

Irrigation 

10 

79 

Chamaraj asagar 

82 . OC 

N&I supply 

11 

81 

Manchanab el e 

39,00 

Irrigation 

12 

83 

Byramangala 

17.00 

Irrigation 



Pig. 5.4 Existing project; 



SALm 
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1 


Pro 1 ect 
Kittur an i cut 


Details of Pig. 5.4 
Existing Projects 

I'ain River 


2 Kudlur anicut 

3 Haluvagilu anicut 

4 Chandravalli anicut 

5 iodihalli anicut 

Hemavathy 

6 Shan]-:artheertha anicut 
7&3 Sreeramadevaru anicut 
9&1B^ 'Mandagere anicut 

10 liemagiri anicut 

11 Kallahalli anicut 

12 Chiklihole anicut Chilclihole 

13£:23^ Ramanathapura anicut 
14 Kattepur anicut 

Cauvery u/s of 

15c.:17 Charaaraja anicut Krishnarajasagara 


16 Mirle anicut 


18 RamasamudraiTi anicut 

19 Tipp.ur anicut 

20 Hanagode series anicut 

21 Katte raalslawadi anicut 

3 uv a r n a va thy 

22&23 Siriyur canal anicut 

24 Marachalli Awandur 
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Detail 

s of Fig, 5. (...Contd 

25 

Devaraja anicut 

2 S , 27 
& 4E^ 

Viriginadi anicut 

20 

Yadathittu anicut 

29 

Bangaradoddi anicut 

3 0 

Mahadevapura anicut 

31 

Ramasv/aray anicut 

32 

Raj apararaeswari anicut 

33 

Madhavamantri anicut 

34 

Lalcshmanapiara anicut 

35£:36 

Kullahalli anicut 

37 

Konglawadi anicut 

3 8&39 

Saragur anicut 

40 

Alur anicut 

41 

Mura.lahalli anicut 

42 

Hosahalli anicut 

43 

Bandilcere anicut 

44 

Gundal anicut 

45 

Maddur anicut 

46 

Dhangare anicut 

47 

Seven anicuts 

4S 

Sivasamudra anicut 

3Bi 

I'crishnaraj asagar 

5B 

1 

Thaggehalli anicut 


Gauveri^ d/s of I^ishnarajasa 


Ka-bini 


Suvarnavathi 


Giindal 


Shirasha 


Cauvery 



Details of Fig. 5.4 (...Contd. ) 


6B^ Chandanahalli anicut 

73, Marconahalli 

SB-, Kanva 

93^ Byramangala 

10B.J Chamarajasagar 

132 Nugu 

232 Childcahole 

332 Mangala 

4 B 2 Suvarnaraukhi anicut 

5 3 2 Hebb anal 1 a 

20^ Hemavathy 

5C^ Kabini 

9C-, Ilanchanabele 


Shimsha 


Arlcavathi 


Kabini 


Suvarnavathi 

Shimsha 


trbavatl'ii 

Kabini 

Kemavathy 


ICaDini 

Arlcavathi 
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The first dam constructed v/as Krishnarajasagar in the 
year 1931, In the last tx-io decades, 11 more projects are 
constructed. Obviously no interdependencies were considered 
in planning and designing these projects. Besides, there are 
a number of minor irrigation schemes consisting of diversion 
\-7orlcs and tanks in the basin. The details of these schemes 
are given in fig. 5.4. 

5,10 Potential Projects 

On the basis of India heterological Department (IMD) 
identification of drought areas, Karnataka has nearly 21,500 
sq.j-cm of such areas in the Cauvery basin itself (irrigation 
Commission, 1972). Thus, major portion of the basin area is 
drought affected. Therefore, it is necessary to augment the 
irrigation supplies by constructing projects to impound v/ater 
and exploit the available ground water potential. 

There are tx* 7 enty one potential sites where dams can be 
constructed and also there are 3 locations for well fields. 


In the present analysis, it is assumed that at each site 
three different scales of projects can be constructed and each 
scale of development is treated as a separate project. Table 5. 
summarizes the details of all potential surface water projects 
and table 5.9 gives the details of ground water projects (fig. 



Fig. 5.5 Potential projects 


37 



Details of Fig. 5.5 


Nev7 Projects 


1C. 


3C. 


SC 


'■'^1 

8C^ 

IC 2 

SC. 

4C^ 


IOC2 

IIC2 

I 2 C 2 

13C. 

14C. 

15C. 

16C. 

ISC. 

4 . 

19C. 


Vetehole reservoir 
Haranoi reservoir 
Taralia reservoir 
Sagare docldaliere reservoir 
Suvarnavathi reservoir 
Gunda 1 res er vc i r 
Yagachi reservoir 
liosapatna reservoir 
Chilclihols reservoir 
La’csiimanatheer tha reservo ir 
! '"udr ec/un d ih alia res er vo i r 
Upper Sh-imsha reservoir 
Uallurv. amanil'.ere reservoir 
Iggalur reservoir 
Hebbadahalla reservoir 
Belthur reservoir 
Changavjadi reservoir 
Doddihalla reservoir 
Uduthorehalla reservoir 
Arkavathi reservoir 
Devalapura reservoir 
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Table 5.3 

Potential surface reservoir projects 


Sl.Kc 

) . No de No . 

Project 

Maximum 
capacity 
(IOC Kect-m) 

Ranarks 

1 

1 

Hesapatna 

1 45 . 00 


2 

2 

ChikliholG 

6.00 


3 

5 

Harangi 

183.00 


4 

9 

Lalcshnianathcertha 

19.00' 


5 

11 

Vctehole 

23.00 


6 

13 

Yagachi 

27.00 


7 

30 

Tar aka 

7 4.00 


8 

37 

Sagaradoddalcere 

6.00 


9 

43 

Kudregundiliall a 

8.00 


10 

47 

Devalapura 

7.00 

All are irr 
gation 

11 

48 

Nalluruamanikere 

29.00 

pro j ects 

12 

53 

Suvarnava th i 

31.00 


13 

55 

Hebb ad aha 11a 

4.00 


14 

62 

Upper Shimsha 

57.00 


15 

69 

Iggalur 

6.00 


16 

76 

Gundal 

3 0.00 


17 

S6 

Arkavathi 

42.00 


18 

89 

Belthur 

5.00 


19 

90 

Uduthorehalla 

22.00 


20 

93 

Changawadl 

8.00 


21 

95 

Doddihalla 

4.00 
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Table 5.9 

Ground '-jator potential 


SI. 

ho. 

Node No. 

Location 

(Hydrological unit) 

Max. capacity 
(100 Hect-m) 

1 

16 


90,00 

2 

64 

and 

■226-. 00 

3 

91 

r 

136.00 


5.11 Deterrriinistic Model 

The inflows to the systaa are assumed to be known with 
certainty. Only mean annual flovjs are taken in the analysis/ 
since extrari-e event design for one purpose may not bo suited for 
other purposes. Input to the system at different points are 
talcen from historical records. 

As mentioned earlier at every site different scales of 
developments are talcen as separate projects. There are sixty 
projects in total, i.e. 12 existing projects and forty— eight ne^’ 
projects. 

With economic and physical characteristic of the systoa 
outlined, the objective function to be maximized is; 

Maximize the present value of net benefits, Z 

where 
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3P 12 30 60 


3 0 60 

I \ Z g, 

t^=l s 


-JO -’^^_ JvL. 

=13 s s 


30 12 

I \ I o 

t=l S ]c=l 


3.0 


- I \ 


Of) 


-1 


•> O O 

4, ^-kt kt 


k=13 


(5.11.1) 


where is a function of irrigation water suddIv, cower surnO 

s - - . , . 

and municipal and industrial water supolies, i.e., - 

£(451^^^ + 96 + 1000 t = 1,2 ^ 

1 = Kharif 

2 = Rabi 

Budgetary constraints 
60 

I 1 10,000,000 for all t„ 

k=13 ^^^s ^ 


Institutional constraints 


60 

y gi x. < 1 for all t,-. 

•j “ Ktl — ^ 

k=13 s 

In any year only one project is introduced. 


30 


I 

t=l 



< 1 


s 

where k' = number of different scales of project that can be 
constructed at any site. For example, at node 1, one of the 
projects, 13 & 14- may be constructed (if at all to be construicted) 
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30 

J, "-t i 1 

t=l s 
s 

A project can be constructed only once. 

Contingent projects 

The project at node point 13 is contingent iv’ith the project 
at node point 12. 

The constraint is v/ritten as 

Q ->- < cf , t 

^23/^4 "s- -^21/22 ""s 

Demand constraints 

The demandi for via ter comprises of irrigation demands power 
demand and municipal and industrial water supply demand. Most 
of the projects in the basin are irrigation projects. The 
projects at node points (fig. S.S) at 21 and 35 are multipurpose 
projects supplying v/ater for irrigation, power and municipal and 
industrial supply. The project at node point 79 is only municipal 
and industrial water supply project. Hence dem.and constraints 
can be written as 

X, , > I, , for t = 1,2 

jCl. iCC 

and k for all projects except projects at 
21, 35 and 79 

X, , > I, , + + VD for t = 1,2 

let — kt let kt 

and for k = projects at 21, 35. 

for t = 1 , 2 

and It = 7 9 « 



It can be seen from the fig. 5.3 that some of the agricul- 
tural districts are supplied with water from more than one 
project. These irrigations are jointly constrained. 

^ K = set of ore jects stioplying water to 

kek - s 

s ...... 

some irrigation district 

AL^ = maicimum irrigable land with 

supply from 1:^ projects expressed in 

the unit of v/ater. 

i.e. irrigation districts 6 , 14, IS, 25, 26, 34, 45, 54, 65, 32, 
92 and 96 are supplied with v/ater from more than one project. 

Maximum storage constraint 


S, , < V, , for all ic and t 

ivtl ** ivt 

Minimum storage constraint 

for Ic = 21, 35 and 79 

iCtl. — ‘ iCtl. 

t = t^ and t 2 

Continuity equation 

\t = ^lct-1 + ^ict " ^Ict '' 

and t = t^ & t 2 

Canal constraints 

X - < H for all arcs 

tm — m 


for all arcs 
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The values of and are given in table (B-2) 

Continuity equation at each node. 

At each of the nodes total inflov- into the node and total 
outflovj from the node must be equal 

I - I X = O for i = 1,2,...50C. 

and 

g. = O or i for all jc and t . 

-vL-s s 

5.12 Stochastic Model 

As mentioned in Chapter 3 only natural inflov/s to the 
system are taken as random variables. The probability density 
coefficients for all stream are obtained using Kim's (1968) 
method. This method is used to describe the flow level proba- 
bility. The method consists of deriving from annual stream 
flow data discrete points. From the historical data mean 
and standard deviations are calculated. The minimum annual 
flow is taken as first discrete point. The sixth discrete 
point is the maxim-um annual stream f lovr. The remaining points 
are obtained by ad.ding to the preceding point, the quotient of 
the difference of maximum and minimum annual flov/s divided by 
five . 

Then the probability density coefficient for each interval 
is obtained by the following expression 
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R. - R R. - R 

Probability density coefficient = F ( ^ - F ) 

o c: 

“"R R 

cor i =1, 2/ 3 6 

Rj^ = discrete point 
R = annual mean stream' flov; 

= -Standard deviation 

and F = normal cumulative distribution function. 

From these probabilitA^ densities, cumulative distribution 
function is obtained. 

The cumulative distribution graph is obtained for all 
streams. Since this involves a large amount of data, only a 
few representative cumulative distributions are given (fig. 5.6 

In the analysis 75 percent probability is used for the 
constraints involving stochastic variables. 

The details of chance constrained formulation with linea^r 
decision rule is given in Chapter 3, In the present chapter 
the application of this formulation is briefly presented. 

As mentioned before, reservoirs at node .points 1, 2, 5, 9, 
13, 30, 35, 37, 42, 43, 47, 49, 53, 55, 53, 62, 56, 71, 76, 79, 
89, 90 and 95 are independent reservoirs since inflows to these 
reservoirs are the only natural flows. Let these reservoirs be 
denoted by a siib-set G 

i.e., G = (1, 2, 5, 9, 13, 30, 35, 37, 42, 43, 47, 49, 53, 55, 
53, 62, 66, 71, 76, 79, 33, 89, 90, 95). 




FIG.5-6 CUMULATIVE DISTRIBUTION CURVE 
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The inflovTs to the reservoirs at nodes 17, 21, 63, 69, 31, 
86 and 93 consist of natural flo".7S plus certain portion of 
releases from the respective upstream reservoirs. These 
reservoirs are called dependent reservoirs and denoted by a sub- 
set H, 

i.e., H = (17, 21, S3, 69, 31, SC, 93). 

Kence the complete set of reservoirs is the union of G and H 
sub-sets. 


The linear decision rule for the independent reservoir 
sifo-set G will become (Refer equations 3.3.6 and 3.3.8), 


Gt 


'Gt 


Gt 


for G (1,2,5.. . ,95) 


and t = t^ & t 2 


X 


Gt 


^Gt-1 ■'’"Gt-1 '^Gt 


Similarly for dependent reservoir suo-set K the linear 


rule 

(Refer 

equations 3.3.12 and 3.3.18). 


^Ht 

" ^Ht '"Ht 

and 

^t 

^ ^It + ^Ht-1 “ ^""Ht J,, C%t^^'Gt-l+^ 

"'"up 


for all H (17, 21, ... 93) 


decision 


Gt-1 



t = t^ & t2 


where = respective u/s reservoirs. 


Equations (3.3.13 & 3.3.16) are used to convert stochastic 
constraints to equivalent deterministic constraints. 



The minimxjm release constraint takes the form for indepen- 
dent reservoir set G 


^ ^^Gt ^Gt-1 ^ 


Gt: 


^°Gt “ -Wl + -V-t^ 1 (1 - 0.75) 


D - 4- b^, < F' 

wti, i_ — 1 'Of in — . 


(0.25) 


for dependent reservoirs sub-set 


°Ht + <''kt bt - fkt-l b.t-l’ ^■bo.25) 


up 

Similarly minimum storage constraints for independent 
reservoir sub-set G vrill be 


jnin 

^Gt 


b_ < F ^ (0.25) 

vj u, — K 


for dependent reservoir siib-set H 


^Ht %t - ^R 

The normal distribution is talcen for all reservoir inflows- 


5.13 Model with Evaporation Losses 

Since the basin lies in sani-arid region the evaporation 
losses from water surfaces are quite considerable. As mentioned 
in Chapter 3 the evaporation losses are given by an expression 
- 6^ (S^,j. -i- S, , T ). These losses are estimated using 

let Jet jC u 

meteorological data, since available data is too meagre. The 
area-capacity curves are obtained using elevation-area and 
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elevation— storage curves v/hicli are available in project 
reports of different dams (fig. 5.7a,5.7b), The evaporation 
parameter foj; Kharif and Rabi seasons is obtained by 

fitting a straight line to the losses at different storage 
levels. The values of are given in table 5,10. 


Table 5.10 

Values of evaporation parameter 


SI. 
Ho . 

Node No = 
(Fig. 5,5) 


Values of 


Reservoir 

Kharif 

season 

Rabi 

season 

1 

2 

3 

4 

5 

1 

1 

1 

Hosapatna 

0,030 

0.100 

2 

2 

Child ihole 

0.030 

0.100 

O 

5 ■ 

Karangi 

0,035 

0.120 

4 

9 

Lakshmanatheertha 

0.035 

0.120 

5 

11 

VeteholG 

0.035 

0.120 

6 

13 

Yagachi 

0. 040 

0.125 

7 

17 

Hern ava thy 

0. 040 

0.130 

8 

21 

I'-jrishnar a j asagar 

0. 045 

0.130 

9 

30 

Tara'-ia 

0. 040 

0.120 

10 

31 

Hebballa 

0.040 

0.110 

11 

35 

Kabini 

0.040 

0. 100 

12 

37 

Sagaredoddalcere 

0.040 

0.120 

13 

42 

Nugu 

0.045 

0.120 

14 

43 

Kudr egundi hal 1 a 

0.045 

0.125 
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Table 5.10 (...Contd.') 


i 

2 

3 

4 

5 

15 

47 

Devalapura 

0.045 

0.125 

1 6 

49 

Talluruamanilccre 

0.046 

0. 125 

17 

53 

S u var n a va th i 

0.047 

0.125 

18 

55 

Hobb adalial 1 a 

0.047 

0.120 

19 

58 

Childcahole 

0.047 

0.120 

20 

62 

Upper Shimsha 

0.050 

0.130 

21 

63 

Marconahalli 

0.050 

0.130 

22 

66 

Hangala 

0.050 

0.130 

23 

69 

Iggalur 

0.050 

0.125 

24 

71 

Kanva 

0.051 

0.130 

26 

79 

Chamaraj as agar 

0.052 

0.130 

27 

81 

M anchanabel e 

0.052 

0.130 

28 

83 

By ramangal a 

0.051 

0.130 

29 

86 

Arkavathi 

0.052 

0.130 

30 

89 

Belthur 

0.055 

0.140 

31 

90 

Udutho rehall a 

0.055 

0.150 

32 

93 

Changawadi 

0.060 

0.150 

33 

95 

Doddihalla 

0. 060 

0.150 


0 

4 

8 12 16 20 

Volume (100 hect-m) 

24 

28 

32 



F1G-5-7(a) AREA VOLUME 

CURVES 





oc 


6C 


>X 


10 


0 


0 2 • 4 


. 6 8 10 12 14 16 10 

Volume (TOO hect-ml 


FIG.5-7{b) AREA VOLUME CURVES 
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In the stochastic model with evaporation losses the 
objective function and the constraints except the minimum 
release and the minimum storage constraints remain same as 
those of previous model. Hence complete model is not repeated. 
However, the minimum release constraint and minimum storage 
constraint are modified using the values of i 

From section 3.5 the linear decision rule for 

independent reservoir sub-set 3, (G = 1, 2, 5, ...) will be 


X 


Gt 


%t-l 


^ ^"Gt-1 '^Gt-I ^ ^ 


?t 


and S 


Gt 




Then the minimum release constraint v/ill be 




^t 

6_. . '‘^Gt-1 ' ■■'Gt-l 


(R. 


Gt-l 


+ - ^Gt^ - 


e. 


(l^Gt + °Gt> ^ - "-Gt-l 

L- 


:< F ^ ?0.25) 


The minimum storage constraint will be 


p > 0.75 


for t =1,2 

G = (1,2,5. 


p , gmin, , o_,5 

>■ /> 5^^^ — 


e. 


p e .. < 0.25 

^ ^ Gt Gt Gt - 


rain q _ 7^ ^ p- 

Gt Gt Gt — 


_i 


(0.251 
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For dependent reservoir siib-set H (H = 17,21,.,.') the 
linear decision rule v/ill be 


X. 


and 


Ht 


Ht 


HHt 

^t-; 


1 + 


Idt-l 


"Ht-l] 


t I 

Gelc 


’^Gt ^Gt 


Ht 


up 


+ 1 : 


Ht 




■It 


"iith these values for snd the minimum release and 

minimum storage constraints takes the form of 


{ %t ^ 


6 


?ek 


Gt Gt 


J Hr 


Ht-1 


up 


Ht 


^ ^ (0.25) 


and 


„min e < „-l 


.T, - b 1. F~ (0.25) respectively. 


Ht Ht Kt 


5.14 Data Details 

The above models are solved using the methodology explained 
in Chapter 4. Table 5.11 shows the location, physical and econo- 
mic data of projects. Any of the different scales of development 
may be built at each location. It can also be seen that projects 
13 or 14 may be constructed at site nirmber 1 and 15 and 16 at 
site number 2, etc. 

5.15 network Details 

As mentioned earlier, a year is divided into two periods 
Appendix B gives the details of all arcs. 


Kharif and Rabi. 
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Table 5. IX 



Physical 

and economic 

data of projects 


SI. 

No. 

__ 

Node Ko. 

Max . 

capacity 
(ICC Kect-m) 

Investment 
(E;. 106 ) 

Maximum 
annual return 
(P'. 105) 


2 

3 

4 

5 

Exis 

1 

ting projects 
17 

960.00 



2 

21 

1140. 00 



3 

31 

11. 00 



4 

35 

541.00 



5 

42 

122. 00 



6 

50 

10.00 



7 

63 

64.00 



o 

66 

7 . 00 



9 

71 

13.00 



10 

79 

82 . 00 



11 

81 

39.00 



12 

83 

17.00 



I'lew 

projects 




13 

1 

145.00 

227.00 

14.50 

14 

1 

110. 00 

200.00 

11.00 

15 

2 

6. 00 

15.40 

0. 60 

16 

2 

3 . 00 

”! C OC 

0.30 

17 

5 

183.00 

110.00 

13.30 

18 

5 

140.00 

32.00 

14- 00 

19 

9 

19.00 

15.50 

1.90 

20 

9 

12.00 

9.00 

1.20 

21 

11 

23. 00 

19.90 

2 . 30 

22 

11 

20. 00 

13.50 

2.00 

23 

13 

27.00 

71.20 

2.70 

24 

13 . 

20.00 

55.00 

2.00 


106 


Table 5.11 (... Ccntd .) 


1 

2 

3 

4 

5 

25 

16 

35 . 70 

35.60 

6. 00 

26 

16 

30.00 

31.50 

4. OC 

27 

30 

74 . OC 

39.50 

7.40 

28 

3 0 

60.00 

37.00 

S . 00 

29 

37 

6 , 0 C 

9 . 00 

0.60 

30 

37 

3. 00 

8.00 

0.30 

31 

43 

8 . 00 

14.00 

0 . 30 

32 

43 

4.00 

10.00 

0 . 40 

33 

47 

7.00 

16.00 

0.70 

34 

47 

4.00 

12.00 

0.40 

35 

49 

29.00 

19.00 

2.90 

36 

49 

15.00 

12.00 

1.50 

37 

53 

31.00 

24.80 

3.10 

38 

53 

25.00 

16.20 

2. 50 

39 

55 

4.00 

9.30 

0.40 

40 

55 

2.00 

7.00 

0 . 20 

41 

62 

57.00 

62.00 

5.70 

42 

62 

40 . 00 

52.00 

4.00 

43 

64 

20 . 00 

31.00 

6 . OC 

44 

64 

15.00 

18.00 

4.00 

45 

69 

2.00 

18.00 

0 . 20 

46 

69 

1.00 

13 . 20 

C . 10 

47 

76 

30.00 

29.60 

3.00 

48 

76 

20 . 00 

2 2.25 

2 . OC 

49 

36 

42 . 00 

37.65 

4 . 20 

50 

86 

30 . 00 

20.25 

3.00 

51 

89 

5.00 

14.20 

0.50 

52 

89 

3.00 

11.00 

0.30 

53 

90 

22.00 

46.80 

2.20 
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Table 5.11 (. . .Contd.) 


1 

2 

3 

4 

5 

54 

90 

15. CO 

38.35 

1.50 

55 

91 

20.00 

25.00 

3.50 

5 6 

91 

15.00 

15.00 

2.00 

57 

93 

8. 00 

8.60 

O. 80 

53 

93 

5. OO 

6.20 

0.50 

55 

55 

4.00 

7.70 

0. 40 
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upper bound, lower bound, value of return from one unit 
(1 Hect-m) of vrater through each arc of the system and initial 
f lovrs . 

The inflow to the system, (natural stream flows arcs 1 
through 40 for i.^iarif season and arcs 225 through 264 for rabi 
season) and original reservoir levels (arcs 41 through 76) 
originate from source node 299. The flov/s in arcs, i.e, arcs 
1 through 40 arcs 225 through 264 and arcs 41 through 76 are 
assured by setting the upper and lovrer bounds of these arcs at 
the same value. All exit flov;s from the systea (arc 403 through 
485) and all final reservoir voluraes (arcs 377 through 412 are 
removed to the sinlc node 30C. The artificial arc 437 connects 
the sinlc nocte to source node 299 to complete the circulation 
of network. 

Stream flov/ records are availalole for all the streams. 

The length of the records are not the sam.e for all the streams. 
Since it is a volirainous data, listing all the data is not 
feasible. Only annual mean flov/ values and 25 per cent probabi- 
lity values are given in table 5.12. 

5.16 Demands for Water 

The analysis is considered under the basic demand pattern 
of irrigation, power and municipal and industrial supplies. 

The danands for power and municipal and industrial (Table 5.12) 
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Table 5.12 


Stream flow data 
(lOO Hect-m) 


SI. 

Stream 

Kharif 

Rabi 

Ho. 

25 per- 
centile 

50 per- 
centile 

25 per- 
centile 

50 per- 
centile 

1 

2 

3 

4 

5 

6 

1 

Cauvery 

990 

1220 

75 

ICO 

z 

Chil'lihole 

52 

65 

4 

6 

■J 

Harangi 

410 

520 

22 

30 

4 

Lalcshirianatheertha 

67 

90 

'I 

! 

10 

5 

Vetehole 

68 

81 

6 

o 

6 

Yagachi 

320 

405 

21 

25 

7 

Hemavathy 

1130 

1400 

240 

300 

o 

u 

Cauvery 

ICri shnara j asagar 

1440 

1800 

310 

400 

9 

Betv/een Fx ishnar a j asagar 830 
and Kabini 

HOC 

110 

150 

10 

Tar aka 

120 

160 

14 

20 

11 

Hebballa 

68 

80 

3 

/!. 

12 

Kabini 

1500 

190C 

162 

200 

13 

Bandiganda 

(Sagare Doddcilcere) 

33 

4-7 

6 

10 

14 

Nugu 

430 

676 

33 

44 

15 

Kudregundi Hall a 

40 

54 

s 

6 

16 

Devalapura 

23 

32 

2 

3 

17 

G-undlu 

(Nalluru 3>raani Kere) 

80 

115 

7 

10 

18 

S uvar n ava thy 

71 

95 

12 

20 
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Table 5.12 (...Contd.) 


1 

2 

3 

4 

5 

6 

19 

iiebbahall a 

8 

10 

0 

2 

20 

Chikkahole 

81 

124 

2 

6 

21 

Upper Shimsha 

550 

7 40 

60 

SO 

22 

Shimsha 

(Marcona Hally) 

140 

185 

21 

30 

23 

Hangala 

24 

34 

3 

6 

24 

Shims.ha 

(Iggalur) 

300 

430 

5 

10 

25 

ICanva 

30 

440 

3 

5 

26 

Between I'abini & Shimsha 

270 

360 

29 

40 

27 

Gundal 

17 

25 

0 

2 

28 

Arlcavathi 

(Charaarajasagar ) 

64 

85 

6 

10 

29 

Arkavathi 

(Manchanabel e ) 

60 

80 

5 

10 

30 

Vr i shabhava thi 
(Byramangala) 

24 

30 

2 

5 

31 

Arkavathi 

150 

220 

7 

lO 

32 

Shimsha to Arkavathi 

130 

180 

13 

20 

33 

Thattehal 1 a 
(Belthur) 

25 

38 

2 

5 

34 

Uduthorehall a 

46 

60 

6 

lO 

35 

Thattehal la 
(Changawadi) 

15 

20 

1 

3 

36 

Doddihalla 

9 

12 

1 


37 

Arkavathi to Boundary 

120 

150 

21 

3C 
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Table 5.13 

Power and Municipal & Industrial Denand 


SI. 
Mo . 

Use point 

Demand 
Khar if 

(100 Hect-m) 

Rabi 

Source 

, 

Power demand 





1 

28 

1000 

1000 

21 


2 

36 

250 

250 

35 


Msa 

demand 





3 

23 

15 

15 

21 


4 

27 

15 

15 

21 


5 

80 

83 

82 

35,79 



*1#' 
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supplies are taken as mandatory. These demands are met by 
specifying tipper and lower bounds for arcs 181 through 184 
(power) and arcs 185 through 188 (M&I supply) for kharif season; 
and for arcs 369 through 372 (power) and arcs 373 through 376 
(M&I supplies) for rajci season, at the same level. 

The demand for irrigation is asstmed to be increasing at 
a specified schedul.e. It may be decided to bring certain 
hectares of arable area under irrigation under each project per 
year. These demands are introduced in the following statement 

HI (arc) = HI (arc) + KI*(t - t^) 

This statement gives the role of increase of irrigable area KI 
and at what year ti^ the actual irrigation starts after- the 
particular project is introduced into the system depending upon 
management policy and prevailing conditions. For each irrigation 
district one such statement is introduced. The maximum limit is 
specified by the statement of the following type 

IF (HI (arc).GT.IR) HI (arc) = IR 
where the IR is maximum irrigation facility under the project. 

Table 5.13 gives the details of irrigation demands in 
Iharif and irabi seasons of the year. 


113 


Tabl e 5.13 

Irrigation denand 
(lOO Hect-m^ 


SI. 

No. 

Irrigation- 

district 

Season 

IQiarif 

Rabi 

Source 

1 

2 

3 

4 

5 

1 

4 

175 

115 

1, 2 

2 

6 

259 

123 

5,19 

3 

10 

16 

16 

9 

4 

12 

30 

20 

11 

5 

14 

100 

44 

13,16 

6 

18 

260 

140 

17, 16 

7 

20 

280 

193 

8 

8 

22 

650 

535 

21 

9 

25 

1054 

650 

17, 62 

10 

26 

500 

260 

21, 63 

11 

29 

1000 

1000 

24 

12 

32 

16 

9 

31 

13 

3 4 

570 

230 

30,35 

14 

3S 

6 

/t 

37 

15 

45 

120 

72 

44 

16 

46 

ICO 

100 

41 

17 

48 

9 

6 

47 

18 

50 

20 

9 

49 

19 

54' 

360 

235 

35, 53 

20 

56 

o 

2 

55 
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Table 5.13 ( Contd. ) 


1 

2 

3 

4 

5 

21 

59 

12 

S 

CD 

m 

22 

65 

25 

25 

66, 64 

23 

68 

150 

150 

67 

24 

70 

20 

16 

69 

25 

73 

20 

14 

71 

26 

75 

650 

550 

74 

27 

77 

8 

4 

76 

28 

82 

25 

15 

81, 64 

29 

84 

15 

15 

83 

30 

87 

40 

20 

86 

31 

92 

32 

20 

89,90,91 

32 

94 

16 

9 

93 


33 


96 


4 


2 


95 


Chapter 6 


RESULTS, Discussions mD CONCLUSIONS 


As mentioned in Chapter 5, the follov/ing cases are 
considered in the analysis: 

a) deterministic streara flows 

(seasonal mean flows) 

b) stochastic streara flows 

c) model v;ith evaporation losses 

d) model vrith retiarn flows. 

All the above cases are studied for two discount rates 
namely, 3% and 4%. The operation, maintenance and replacement 
costs are tahen into account as a certain percentage of capital 
costs such as 0.25 percent and 0.50 percent. 

The systera is a very’ large one, consisting of 33 surface 
reservoir projects and 3 well field projects. The network (of 
reservoirs, well field units, cane^ls, river channels, irrigation 
districts, municipal and industrial use points and power houses) 
consists of 300 node points and 437 arcs in the case of two 
season model. Because of the sire of the model and number of 
combinations studied, no atterapt is made here to present all 
the resiilts for all the cases considered. The results of 
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c-eterministT c and stochastic models only are given in Appendix B. 

A list of selected projects and their sequence of intro- ' 
duction into the system based on deteirministic flows using a 
discount rate of 3 % and a benefit rate of Rs. 45/hect-m for 
irrigation is given in table 6,1. It can be seen from this 
table that all the reservoirs except upper Shimsha are selected 
to their maximum capacity. The levels of the reservoirs in the 
system at the end of different seasons are presented in table 
6.2. The amomt of water directed to irrigation districts in 
different seasons in case of deterministic and stochastic models 
is given in table 6.3 (see also tables B-4 and B-6). The 
irrigation districts whose recruirements are not met fully are 
underlined. 

It can be seen that in hharif season, in the case of both 
deterministic and stochastic models, the requirem.ents of irriga- 
tion are met except in only one irrigation district (25). But 
in the case of rabi season, requirements are not fully met in 
many irrigation districts. For exarnple, for the deterministic 
case, requirements of 25 are not met and in the stochastic case 
requirements of 25, 26, 29, 32, 34, 54, 59, 68, 70, 73, 84, 92, 
are not met. In the case of stochastic model, the final network 
configuration is given in table B— 6. Prom the flow values in 
different arcs of the deterministic case with 3 % discount rate 
(Table B-4) , It can be seen that there is no probability of 
impounding more water upstream of Krishnarajasagar because of 
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Table 6.1 

Selected projects and their sequence of introduction into the q'-t3r'. 



( deterministic 

case, d 

iscount ra 

te 3?4) 


31. 

ho. 

Project 

Mode 

Mo. 

Reser- 

voir 

No. 

Capacity 

(ICC 

Kect-m') 

Year of 
intro- 
duction 

Motation 

1 

2 

3 

A 

5 

6 

7 

1 

Sagaredoddahere 

37 

29 

6.00 

1 

29, 1 

2 

Well units 

64 

43 

20. 00 

2 

43, 2 

3 

Yagachi 

13 

23 

27 . 00 

3 

23, 3 

4 

Upp e r Sh ira sh a 

52 

42 

40.00 . 

4 

42, 4 

5 

Karangi 

5 

17 

1 83 . 00 

5 

17,5 

6 

Well unit 

16 

dl D 

35.00 

6 

25, 6 

7 

Uduthor shall a 

90 

53 

22.00 

7 

53/7 

3 

Taralca 

3 0 

27 

74.00 

8 

27,3 

9 

Arkavathi 

36 

49 

42. GO 

9 

49,9 

10 

Hosapatna 

1 

13 

1 45 . 00 

10 

13,10 

11 

Gundal 

76 

47 

.30. 00 

11 

47, 11 

12 

Suvarnavathy 

53 

37 

31.00 

12 

37,12 

13 

Well unit 

91 

55 

20. 00 

13 

55,13 

14 

Vetehole 

11 

21 

28.00 

14 

21,14 

15 

Iggalur 

69 

45 

2. 00 

15 

45,15 

16 

Kallur Amanikere 

49 

35 

29.00 

16 

35,16 

i *7 

X ! 

Child ihole 

2 

15 

6.00 

17 

15,17 

13 

Lakshmanatheertha 

9 

19 

19.00 

1 

19,18 
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Table 6.1 (...Ccntd. ) 


1 

o 

3 

4 

5 

6 

7 

19 

Devalap-^.ira 

47 

33 

7. 00 

19 

33, 19 

20 

Kudr egundihal 1 a 

4-3 

31 

3.00 

20 

31,20 

21 

Bel thur 

39 

51 

5.00 

21 

51,21 

22 

Keb’s adah alia 

55 

39 

4. 40 

22 

39, 22 

23 

Changawadi 

93 

57 

s.oc 

23 

57, 23 

24 

Doddihalla 

95 

59 

4. CO 

24 

5 c g 


1 0 

Present value of net benefits = T:;. 0.11206 x 10 

(irrigation benefit rate Rs. 45/Kect-m 
Power vrater rate Rs. 96/'Iect-m 

M&I Water supply rate Rs. 2 000/Kect-m 

*In this column the first fig’ore in each rovr indicates the 
project no. and the second figure indicates the year of 
introduction. For example, 29,1 means, project 29 is 
introduced into the system in the year 1. 
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Table 6.2 



Reservoir 

contents at 

(in 100 ! 

the end of the season 

lect-m) 

Kode No . 

.ivxicir..L. ^ It:- / 

Rabi (t,-,) 

Project 

1 

2 


4 

1 

ICC 

145 

I-'osapatna 

2 

9 

6 

Chikihcle 

5 

90 

133 

Harangi 

9 

' 0 

19 

Lal'shmanatheertha 

11 

12 

28 

Vetehole 

13 

8 

27 

Yagachi 

16 

20 

■ J 

VJell unit 

17 

300 

960 

Hemavathy 

21 

600 

1 040 

Kri sbnara j asagar 

30 

40 

74 

Taralca 

31 

6 

11 

iieiO-O alia 

35 

350 

5 41 

Kahini 

37 

2 

6 

Sagaradocidahore 

42 

60 

122 

Tugr, 

43 

3 

8 

Rucr egund ih alia 

47 

3 

7 

Devalaptira 

49 

o 

29 

Rallure amani I'ore 

53 

.0 

31 

Suvarnavathy 

55 

1 

4 

Hebbadahal la 

58 

5 

1C 

ChikTcahol e 
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Tabl e 6.2 ( . . . Contd. ) 


1 

2 

3 

4 

62 

12 

40 

Upper Shimsha 

63 

30 

64 

Marcona Hally 

64 

30 

60 

Hell anit 

66 

3 

7 

Mangala 

c 9 

2 

2 

Iggalur 

71 

9 

18 

Kanva 

76 

p 

30 

Gundal 

79 

.-1 r * 

82 

Chamara j asagar 

B1 

20 

39 

K anchanabel e 

83 

’ 1 . 

17 

Oyramangala 

36 

1 0 

42 

Arkavathi 

89 


5 

Oelthiir 

90 

3 

22 

Uduthorehal 1 a 

91 

10 

35 

Hell unit 

93 

2 

8 

Cnangav;adi 

95 

4 

.d 

Dc-ddihalla 
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Table 6.3 

Quantity of water supplied to various irrigation 
districts against their demands 
(in ICC Kect-m) 

(Sum indicates supply from more than one project) 


SI. 

Ko. 

Irriga- 

tion 

district 

Kharif season 


Rabi season 

Demand 

Supply 

Demand 

Supply 

Deter- 

ministic 

CS-SS 

Stc chas- 
tic 

case 

Deter- 

ministic 

case 

Stochas- 

tic 

case 

1' 

2 

3 

4 

5 

6 

7 

8 

1 

4 

175 

i53-fl7 

15B-fl7 

115 

1 Oc +1 0 

1 05 +8 

2 

6 

259 

49-f21C 

49+210 

123 

33+90 

33+90 

3 

10 

16 

15 

15 

16 

1 6 

16 

4 

12 

3C 

.r'. 

0 V. ' 

30 

20 

20 

20 

5 

14 

ICO 

100 

100 

44 

44 

44 

6 

18 

260 

2 25 -f-3 1) 

225+35 

140 

115+25 

115+25 

7 

20 

230 

280 

28C 

223 

193 

97 

3 

22 

650 

65 0 

65 0 

535 

510 

20 

9 

25 

1054 

626+134 

35 6+3 4 

650 

6004-23 

397-^23 

10 

26 

5 00 

440+60 

25+60 

260 

1 1 3+40 

0+40 

11 

29 

1000 

1000 

1000 

1000 

406 

335 

12 

32 

16 

16 

16 

9 

9 

8 

13 

34 

570 

500-:-70 

500+70 

230 

101+30 

63y-3C 

14 

38 

6 

6 

6 


4 

4 

^ 5 

45 

120 

120 

120 

72 

72 

72 

16 

46 

100 

100 

100 

100 

100 

100 

17 

48 

9 

9 

9 

6 

6 

6 

-t 

1 ,0 

5 0 

20 

20 

20 

9 

9 

9 
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^ a.Dl0 

6.3 

( . - . Contd. 

) 





1 

2 

3 

4 

c; 

w' 

6 

7 

8 

19 

5 4 

360 

300+60 

300+43 

235 

0+35 

0+35 

20 

56 

3 


3 

2 

2 

2 

21 

59 

1 2 

12 

12 

S 

O 

1 . 

22 

65 

25 

4_ ' 

5+20 

25 

5+20 

5+20 

23 

68 

1 c r"\ 
j. .J 

1 K p- 

150 

15 0 

lie 

82 

24 

70 

20 

20 

20 

16 

16 

C 

25 

73 

20 

20 

20 

14 

14 

12 

26 

75 

65 0 

550 

S50 

550 

550 

550 

27 

77 

O 

3 

1 

4 

4 


28 

82 

25 

15 -^-10 

15+lG 

15 

5+lC 

5 +10 

29 

84 

15 

1 5 

15 

15 

14 

14 

30 

87 

40 

40 

40 

20 

20 

20 

31 

92 

57 

32+20+5 

54.20+32 

45 

21+4+20 

4+20+20 

32 

94 

15 

16 

13 

9 

9 

O' 

33 

95 

4 

4 . 

4 

2 

2 

2 
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the requirements of the reaches dovm stream of the dam. It 
implies that water resources are fully utilized. It is the same 
in case of stochastic flows also. 

The valuviS of decision variables, (section 3.4.1) are 
given in table 6. 4. These values give the quantity of water 
that flows Gov.T'. thr? strear'i ah£.nnel3. In rabi season the value 
of for arcs 266, 270, 272, 277, 279, 281, 289, 294, 297, 304, 

3 05 and 318 is zero indicating that these reaches go dry in rabi 
season, i.e., no water flows dov^n the river channels from 
reservoirs Chiklihole, Yagachi, Lakshmanatheertha, ICrishnaraia- 
sagsx, Hebballa, ICaIjini, Devalapura, Suvarnavathi, Cnikkahole, 
Iggalur and Changawadi. It can be seen from the table that the 
values of s for lurishnaraj asagar, Hemavathy, Yabini, S^J.varnava 
Kanva, Gundal, Chamaraj asagar, Byramangala and Uduthorehalla arc 
less than the value of 0.40. Hence, it can be said that the 
available v/ater resoiirces are used efficiently in these cases. 

Dual Variables 

The values of for original network are zero for each arc 
{ = O') and for final network the values are cjiven in table 6.4 

From these values it is evident that the valuesof dual variables 
y. . for 150 arcs (table B- , ) are not aaual to the value of 
corresponding ijnit rate of return per hect-m. 

For original network conf igiiration the value of dual 
variable y. . is equal to unit rate of return, i.e,, from 
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Table 6.4 

Fracuicr.s cf release (H*) %-7hich will flov; down 

the stream 


Project 
at node 


Annual mean 

i’ :'C. ■ G. 1 '* 

flows 

sole) 

R-lbi 

Stochastic 
(25 percen 
rCnarif 

flows 

tile) 

Rabi 

1 


2 


4 

5 

1 

X 


0.335 

0. 275 

C. 82 

0.125 

2 


0.7 20 

0. 000 

0. 660 

0. 000 

5 


r. S' 5 

0.930 

G.S50 

0. 710 

9 


■W * w 

0.153 

0.725 

0. 000 

11 


0.5 40 

0. 2CC 

0. 420 

0. 30 

13 


0.7 40 

C. OCC 

0.670 

0.000 

17 


0 , 1 9 S 

0.170 

0. 210 

0.350 

21 


0.13C 

0. 000 

0.000 

0 . coo 

30 


0.445 

0. 445 

0. 1 86 

0.360 

31 


0.790 

0. 400 

0.350 

0. OC'O 

35 


0.360 

0. OCC 

0. 200 

0 . 000 

37 


0. 830 

0. 670 

0,730 

0. 600 

42 


1 . OGO 

1 . COC 

1 , OOC 

1. 000 

43 


1 . COC 

1 . OOC 

1 . 000 

1 . COO 

47 


0. 670 

0.1 40 

C.5 20 

0.000 

49 


0.730 

0.7 00 

0. 670 

0. 67 G 

53 


0.166 

C.186 

0. OOC 

0.000 

55 


C . 5 7 0 

0. 600 

C. 4-00 

0.330 

58 


O. 900 

0.270 , 

0, 840 

0.000 
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TalDle 6.4 (...Conta.) 


1 

2 

3 

4 

5 ■ 

62 

0.95 0 

0.7 SO 

0,940 

0.740 

63 

C.93C 

0,720 

0.900 

0.67 0 

56 

0* Hi C 

0.450 

0.750 

0. 235 

69 

0. 9 SC 

0 . ooc 

0,970 

0 . 

/ j. 

C. 355 

C. 000 

0.050 

0 . ooc 

76 

X -L 

0.730 

0. OCC 

0, 75Q 

79 

0.1 GO 

c. 120 

0.110 

0. 023 

81 

0.77 0 

0. 35 C 

0.630 

0.670 

83 

0.235 

c. oco 

0. OCC 

0. 000 

36 

0, 340 

0.7 20 

0.720 

0.590 

39 

0. 840 

0,500 

0, 770 

0.200 

90 

0.305 

0,120 

0. OCC 

0. 000 

93 

0. 610 

0.310 

0. 000 

0.000 

95 

0. 60C 

0. 60C' 

0.430 

0. 330 
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section 4.3 for those arcs for which the unit rate of return 
is Ps,45, we car. '»-nrite. 


^ij 

1 

1 

j 

-45 = -45 



= rns-x 

(c. 

-■-45 ) 

m 

z . . 

ij 

r= max 

(C, 

-45 ) 



Thus for these arcs in original nefc-jork system 

y . - = 45 

ij 

For final network system the values of <>• for 150 arcs 
(Table B-5') equal to 45. The value of dual variable is 
calculated for one of these arcs and its significance is 
discussed; 

Consider the arc number 185, 

from table 3-3, o- = 45 and b. . = ICOO 

j-J 

Value of C.^ = 45-O-lOQO = -955 


and 

^'ij 

= max 

(G, 

— C . . J 
IJ 



cr max 

(0, 

955) 


^ii 

= 955 



and 


™ max 

(0, 

-955) 


c 


Therefore, increase in supply in this arc by one -unit 
increases the revenue by R3.955. This can be seen by examining 
the network, configuration aroxmd this arc. The arc 147 compete 
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supply) tor water with, arc 185, Hence increase in 
flow in arc 185 by one unit decreases the flow in arc 147 by one: 
unit. Hence the net increase in revenue is (1000-45) = rs.955. 

Similarly significance of dual variables for other •: 

arcs can be explained. 

The final network configuration indicates the flow in 
the arcs 13 4 & 322, i.e., exit flows from the basin for ICharif 
and rabi seasons are 5090 & 629 xinits respectively, for 
deterministic flow model, (unit = 100 Hect-m) . The corresponding: 
values in stochastic model are 2378 mits and 556 units. Using i 
the stochastic values as criterion for decision, it can be said I 
that the maximum value of mandatory releases without reducing i 
the irrigation potential of the basin are 2378 units in kharif ' 
season and 555 units in rabi season. Any constraint to increase ■ 
these values curtails the irrigation potential of the basin. | 

The differences betvjeen initial netvjork. configuration : 

and final networlc configuration indicate (tables B -2 & B- 3 ) the | 
changes in reservoir conditions and irrigation supply- The : 

values of higher and lower bounds change from zero to the values ; 
shown in the table indicating that the projects are introduced i 
into the system. 

The contingent project constraint is included by imposing 
a condition tiiat Yagachi project is introduced bsfors Vetehole 
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pro j ect is introd'acea into the systsn. Only the sequence of 
introduction of projects into the systern is altered whereas, 
the size of projects remain the same as seen in table 6.5. 

Sensitivity Analysis 

As mentioned earlier several cases are studied. Results 
for each case run into about 30 pages of computer sheets and 
hence it is not possible to include all the results here. 
However, the optimal scheduling and sequencing of projects 
(not including the operation policy) for different cases is 
given in table 6.5. 

It is observed that there is no change in selection of 
projects. This is because, the size of the projects remaining 
same even trough there is a change in sequencing of selected 
pro j acts. 

Conclusions 

It can be concluded that mathanatical programming tech- 
nique combining the branch and bound algorithm ana out-of-Kilter 
algorithm is a pcv/erful tool in the analysis cf optimal 
investment policy and operating policy problems of water 
resources systems- The approach incorporates institutional 
and budgetary constraints such that it provides guidance to 
the policy mal:ers. 
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Table 6.5 (...ConTd.) 


1 

2 

3 

4 

5 

51, 21 

29, 22 

31, 22 

51, 21 

51, 21 


39,23 

39, 23 

39, 22 

"> O 9 0 

57,23 

57, 24 

57, 24 

57, 23 

57,23 

59, 24 

59, 25 

59, 25 

59, 24 

3 9 / Z '-'t 


Present value of 
:iet benefits 
in E3. 


0. 31972 
10^’ 


T. 

X 


,2201 

^,10 


0. 1C2C4 


X IG 


>10 


0.11023 


1C 


.10 


0,31 

X 10 


Xn "til is Ccis 0 iOiroj 0 ct 3t noXs 13 

project at node 11 (Vetahole) 


(Yagachi)is 


contingent with 


the 


m CA 
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The mocel presented in this study can be altered 
suitably tc suit the special conditicns of the region and can 
be employed to ary complex v/ater resources system - mult i- 
purpcse-mul tipro j ec t system. 

In the present studV/ the combinatorial programming 
technique is used for solving the complex water resources system 
taking stochastic nature of inflows and evaporation losses into 
consideration. The chance-constraint programming approach used 
in the model helps to decide the operational policy along i/ith 
optimal scheduling ana sequencing policy, for Icnown levels of 
reliability of service conditions. 

The application to Cauvery river basin is done aaore frera 
a point of demonstrating the application of combinatorial 
programming technique with stochastic nature of inflows and 
evaporation losses taken into account rather than making any 
specific recommendations for the de-ralopment of study area. Ih 
economic parameters and inflow-/ data used for some of the stream : 
in the basin are rather crude. Further, in some cases, the 
data had to be estimated because of pa-Jcity of data and secrecy 
attached to the available data duo to intar-state disputes 
regarding sharing of waters- If complete data is available for 
any basin, then the model can be used for making managerial 
decisions. 
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Suggestions for Future Wcrlc 

Reliable values of objective function coefficients are 
required for the analysis, especially information is needed 
regarding the benefits attached to varicus uses of water and. 
various levels of v.7atcr use. 

The- derrtand pattern has to be established suitably taking 
other related topics such as agricultural, socio-economic 
conditions into consideration. The model should be solved 
taking a month as a season to study the dynamic nature of 
supply and demands. 

Mere vic'crous treatment of the problem, is to consider 
the change of costs and benefits of projects and discount rate 


with time. 



APPEKDIX A 


W.ATSR RATES FOR IRRiaATICE;. XUMICIPAL, 
lAAUSTRIAL & POmR SUPPLIES 


Irrigation Ratos 

The water rar.es arc estimatvod in the follo'wing way : 

liia arcant of consumptive vse is worked out for each 
Crop by Christianson formula. Then talcing farm, efficiency 
(0.7) and ccnveyanca efficiency (0.7) into consideration^ tlv: 
total amount of diver-ion is calculated at the head works. 
These arc given in Table A. i along with the land revenues 
suggested by Irrigation Commission, 1972. 

Table A.l 

Tonount of diverted water 


SI . No . 

Crop 

Amount of water in m 

R a t G> Rs/Hg c t ■ ■ 

1 

Sugarcane 

3.90 

150. OC 

2 

Rice 

2.15 

110.00 

3 

Other crops 

0.S5 

40. 00 
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3y Icnowir.g the extent of areas under different crops 
irrigated by existing projects the following aroa-wise distri- 
eution is ass'jreed: I<ice 50% of the area; sucarcane 30?4 of the 
area and orher crops 20% of tiie area. 


then the averacG water 


cstiitated as follows: 


The rate par r.ect-r for each crop: 


Sugar 



Rico 



0 'bn c IT c iro d s 



51.00 

47 . 20 


The weighted average rata, therefore, is 

38.50 X C.30 - 51.00 x C.5C 47. 2C x C. 20 = 46.50 

Based on the 'a}ocve value, in the analysis two rates, 
namely, Rs. 45/Tect-m. and r%50/Hect-ra are used. 


Municipal and Ix'dustrial Supply Rates 

Municipal and industrial supply rates are worked out usin 
the data of Bangalore Hater Supply & Sev;erage Board (E'"''SSE), 

The supply pattern is divided into a nuraber of types depending 
upon the type of consumer (say domestic, hotels, etc. ) and each 
type is side -divided into different groups depending on amount 
of consumption and suitable weightage is given for calculating 
the rates. The follcwing gives the details of -'■rater rate 


calculation 



Table 


2 
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Type 

1 

supply 

( Regul ar 

domestic supply) 



SI. To 

» 

G.0psriCt 

Group 
ing upon 

Rate 

T-’h fts 

VTeigntag 

■e 7'eighted rate 



ar’oun':: 

# • . 

\ uni 1 

of consume tion 
- ICOC litres’) 



1 



2 

3 

4 

5 

(3 X 4) 

1 



Q:'_ 9 

0.4C 

40?4 

0.16 

2 



25-5C 

0. 10 

40:4 

0.12 

3 



5 0-75 

0.40 

10% 

0. 04 : 

4 




0. 6C 

5% 

0.03 

5 



ICC 

above 0. 30 

5% 

0,04-.. - : 







0.39/10CC ; 
litres 

Type 

2 

supply 

(Gouestic 

supply where no property 

taic is collected) 

1 



C-25 

0, 25 

40% 

0.10 : 

2 



9 9 ». 1 O 

0.3C 

40% 

0.12 ; 

3 



50-7 5 

C. 40 

15% 

O 

o 

4 



75-lCO 

O, 60 

3% 

0.013 

5 



ICC and 

atove C.^0 

2% 

0.016 : 







0. 314/lCCC 
litres : 

Type 

3 

supply 

(corranerci 

.al and other such 

establishments) 

1 



r** ' 

Cu 80 

40:, 1 

0.32 

2 



1 O— 9 O 

1 , OC 

30% 

0. 30 

3 



20-4C 

1.50 

20% 

0.30 

4 



40-100 

2 . 00 

5 % 

0.10 

5 



ICO and 

Qjio've 2.50 

5/0 

0.125 


T J ET ■/ 


1 
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T 3^D 1 S *'4.— 3 

Water rates for Ftower Generation 


T^^pe 

Rate Cr.5. /Hect-m) 

Weightage 
in percent 

Weighted rate 

A, C,D,E 

55.00 

20 

11. ce 

B 

100.00 

20 

20. 00 

P 

f3 . on 

20 

13.30 

G 

45.00 

20 

9.20 

1 

235.00 

10 

23.50 

2 

63.00 

5 

3.40 

3 

290. GC 

5 

■ 14.50 




96.40 


The value of 9S/Kect-m is used in the analysis* 



APFErOIX 


i'able B-1 

details of !'et.v;orl; 


Arcs 

rtemar^is 

I'liarif season 


1 - -l-C 

Jr.flc^n 

41-76 

Reservoir initial conditions 

77-134 

Regulated flox-rs 

135-ia: 

Irrigation flows 

isi-iza 

Power and Il&I flows 

Rabi season 

189-224 

Reservoir initial conditions 

225-264 

Inf ].cws 

265-322 

Kcvsilatecl 

323-363 

Irrigation f lov^s 

369-376 

Fo'-7Sr i- M£:I flov/s 

377-412 

Final reservoir conditions 

413-437 

Si^cit rloviS 
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Table 3-2 
Operating Policy 





Iriginal networ': c 

cnf igior 

aticn 





(..711 ur.ios are 

in ICC' 

Hsct-m) 


Arc 

Source 

Sin^: 

ler.efit 

ret® 

(rs/r:ecz.r:.L-"--. 

beterr.i 

Lower 

bonrib 

nistic 

Plow 

rno de 1 

Ronarl's 

.1 

2 



6 



1 

299 

1 

0 1220 

1220 

12 20 


2 

299 

3 

0 60 

55 

65 


3 

299 

- 

r- 

W* iiC ^ 

5 20 

5 20 



299 

o 

C 90 

9C 

90 


5 

299 

1 1 

0 31 

31 

SI 


6 

299 

13 

C 405 

405 

405 


7 

299 

16 

C 35 

3 5 

35 


8 

299 

17 

0 1400 


1 400 


9 

299 

21 

- 

V.,-' 

1300 

1300 


1C 

299 

2 ■! 

0 1 1 GO 

7 1 pn 

HOC 


11 

299 

30 

C 160 

160 

1 60 


12 

299 

31 

r. po 

30 

30 


13 

299 

35 

0 1 9 0? 

1900 

1900 


14 

299 

37 

C rp 

47 

47 

Inflows 

15 

299 

^ 0 
nr i* 

0 676 

676 

67 5 


16 

299 

43 

0 54 

5 4 

5 4 


17 

299 

47 

0 32 

32 

3 2 


13 

299 

49 

C 115 

115 

115 


19 

299 

53 

C 95 

95 

95 


2G 

299 

DD 

0 11 

1 1 

11 


21 

299 

5 3 

C 124 

124 

124 


22 

299 

62 

C 7 40 

7 /\n 

740 


23 

299 

o ^ 

C 185 

IBS 

185 


24 

299 

64 

0 40 

pu 

40 


25 

299 

66 

0 3 4 

34 

n .1 

-x 


2 6 

299 

69 

0 430 

,4 o n 

430 
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Table 3^-2 (. . .Contd. ) 


1 

2 

0 

4 

5 

6 

7 3 ' 

27 

29 S 

71 

. 0 

4C 

40 

40 

25 

2 9 9 

7 4 


36C 

3 sc 

350 

29 

299 

76 

r- 

■9 c 

25 

25 

0 r'= 

299 

70 

0 

85 

cq 

a..' 

85 

31 

2S9 

Ql 

_L. 

0 

80 

SC 

80 

0 0 

29S 

33 

0 

3 0 

3C 

3 G 

0 ^ 

299 

56 

c 

2 20 

220 

220 

34 

299 

oS 

0 

1 0 

13C 

ISO 

35 

299 

7-4 

c 

33 

33 

38 

36 

299 

90 

V-.' 

6 0 

SC 

6C 

37 

299 

91 

0 

2 0 

20 

20 ^ 

Inf lows 

3 8 

299 

93 

0 

Z 

2G 

20 

39 

299 

95 


12 

12 

12 

40 

299 

99 

0 

i D V. 

1 0 ; 0 

150 

41 

299 

1 

0 

f) 

c 

0 

42 

299 

2 

0 

0 


0 

43 

299 

5 

c 

n 

0 

0 

44 

299 

9 

f", 

0 

c 

c 

45 

299 

11 

0 

c 

c 

0 

46 

299 

13 

0 

0 

0 

G 

47 

299 

16 

0 

0 

c 

C 

43 

299 

17 

0 

3 00 

300 

3 CO 

49 

299 

21 

c 

SOG 

600 

600 

50 

299 

30 

0 

u 

0 

0 

51 

299 

31 

c 

S 

6 

6 

52, 

299 

35 

c 

350 

35 c 

350 

53 

299 

37 

0 

0 

0 

C 

54 

299 

42 

0 

60 

60 

60 

53 

299 

43 

0 

C 

0 

0 

56 

299 

47 

0 

C 

c 

0 
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Table 8-2 (• • .^Contd. ) 


1 

2 

3 


5 

6 

7 

a 

57 

Z i;* 

45 

0 

0 

c 

c 


!D 

299 

5 3 

o 

c 

0 

0 



299 

6 5 

c 

0 

0 

0 


./r 

O V,.* 

23 9 

D Ci 


5 

c. 

5 


5 

25 5 

62 

0 

0 

0 

.r*. 


6 2 

29 5 

c z 

n 

30 

3C 

30 


6.3 

259 

64 

rs 

O 

0 

C 

r.' 


64 

259 

z f, 

f 

3 

85. 

3 


65 

2 o o 

f^G 

C 

0 

c 

G 


j o 

0 c 

i j;-- 

/ i 


9 

9 

o 

Ini tie.]-. 








reservoir 

;;:7 

o n~ .i"'* 

/ o 

f) 

0 

f' 

c 

condition; 

^ P 

299 

"9 

G 

45 

45 

45 

(TO'arif ) 

69 

29 9 

ai 

w 

20 

20 

20 


70 

255 

S3 

(--■ 

G 

8 

P, 


7 1 

25 5 

9 

r-‘. 

V-' 

c 

0 

r" 

V„..' 


72 

29S 

PQ 


0 

c 

0 


73 

252 

93 

(3 

c 

0 

c 


74 

295 

O 7 

0 

0 

0 

0 


75 

295 

3' 5 

0 

0 

. 0 

c 


76 

299 

o - 


r.. 

C’ 

c 


77 

1 

o 

■U 

o 

2 00C 

c 

1138 


7 8 

2 

3 

c 

1 

r\ 

V..' 

63 


79 


(-J 

/ 

c 

2100 

n 

1251 

Regulate.d 

CO 

O 

.5 


w 

87 0 

0 

495 

flows 

81 

/ 

Z 

0 

2500 

t 

V 

1746 


82 

Q 

fy 

r. 

IGO 

0 

30 


33 

11 

15 

C 

ICC 

0 

81 


84 

13 

15 . 

0 

510 


3 37 


85 

15 

17 

o 

650 

0 

468 


86 

17 

19 

0 

2000 


15 68 





O' o 
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Table z~2 (...Tcnrcl.) 


1 

2 

3 



6 

7 

o 

117 

71 

72 

0 

35 

c 

16 


113 

7 2 

7 4 


2200 

0 

1 C.C7 


119 

7 4 

7 3 

P 

9 030 

p 

7507 


12C 

7 "S 

7 5 

c 

■-irD 

0 

21 


121 

7 3 

83 

0 

9225 

0 

7" 27 


122 

79 

81 


O’- 

0 

33 


123 

C= J- 



155 

p' 

SO 


124 

O O 

Obt 

85 


25 

0 

13 


125 

05 


c 

ISO 

c 

93 

Regulated 

flows 

126 

35 

83 

p 

53 0 

r\ 

277 

(Kharif ) 

127 

oo 

99 

o 

1C25C 

0 

7964 


128 

39 

93 

r. 

€C 

n 

30 


129 

9 0 

9S 

r% - 

50 

0 

52 


130 

93 

97 

p 

65 

c 

47 


131 

95 

r\ •~7 

'Z) / 

C' 

25 

c 

11 


132 

97 

3 3 

pi 

1 05 

p 

53 


133 

9 3 

oq 


21C 

c 

110 


134 

99 

IOC 

C 

12430 

0 

8209 


135 

1 

4 

45 

0 

0 

0 


136 

2 

4 

45 

0 

0 

G 


137 

5 

6 

i* rr 

c 

0 

0 


138 

9 

10 

45 

0 ' 

p 

0 


139 

11 

12 

45 

c 

c 

c 

Irrioatic 

140 

13 

14 

45 

C- 

p, 

c 

flows 

141 ■ 

16 

14 

1 ” 

'*X-' 

c 

0 . 

0 


142 

16 

13 


c 

0 

0 


143 

17 

13 

45 

260 

260 

2 30 


144 

17 

25 

45 

C 

V.) 

0 


145 

19 

6 

45 

210 

210 

210 


146 

P 

20 

45 

280 

280 

2 SC 


147 

21 

22 

45 

650 

650 

650 
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Table 3-2 ( . . . Contd. ) 


1 

2 

3 

4 


6 

7 

.r. ^ 

178 

91 

96 

45 

w 

0 

0 


179 

93 

94 

45 

C': 

C: 

0 

Irri^atic 

180 

95 

o 

45 


c 

c 


131 

21 

23 

r 

ICCC 

! GCC 

1 coc 


132 

23 

7 4 

95 

1 coc 

ICCC 

iCOC 

Povrer 

133 

3 5 

36 

f'- 

250 

250 

25C 

Sur.pli/ 

134 

3 5 

39 

96 

25C 

250 

- 2.50, . 


1 ^“'5 

O T, 

Z 

23 

1000 

15 

15 

15 


105 

21 

27 

1 coo 

15 

15 

15 


187 

3 5 

3C 

iccc 

4-C 

40 

4C 

supply ; 

133 

79 

SC 

1 nrn 

J, • w 'v. W 

43 

43 

43 


139 

1 

151 

f. 

0 

G 

0 


190 

0 

15 2 

/-■ 

u 

c 

C 

C 


191 

5 

155 

n 

0 

C 



192 

a 

159 

0 

c 

0 

O' 

i 

193 

1 ' 

161 

0 

c 

c 

0 


194 

i J 

163 

C 

r- 

O' 

G 


195 

16 

1 S 6 

C' 

0 

c 

0 


196 

17 

167 

5 

9 6C 

9 60-' 

960 


197 

21 

17 1 

C 

1140 

1 1 40 

1 1 40 

Initial i 

193 

3 0 

130 

f" 

C 

C, ' 

0 

reserve i 
conditio 

199 

31 

1 3 " 

!_,■ 

11 

11 

11 

I 

2 OC 

35 

135 

f) 

541 

5 41 

5 41 


201 

37 

137 

C 

c 


V..'' 


2 02 

42 

1 92 

c 

122 

122 

122 


203 

43 

193 

'v. 

v-' 

C- 

Q 


204 

47 

197 

0; 


C-' 

C 


205 

49 

199 

0 

n 

w 

0 

.r\ 

V,,.' 


206 

5 3 

2C3 

c 

0 


■ c 


207 

5 5 

205 

C' 

0 

n. 

, 0 
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Table b- 2 (..-Oontd 


1 

4- 

3 

4 

5 

6 

' 7 

8 

2C'S 

55 ■ 

2CS 

0 

10 

10 

10 


2G9 

62 

212 




0 


210 

63 

213 

c 

64 

64 

54 


211 

64 

214 

r 

K. 


0 

G 


212 

66 
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Operating policy 

CDeterTni!lstL''fiTwr 

tlows, discount rate 3%) 


1 

299 

1 

40 

299 

99 

41 

42 

299 

299 

1 

2 

43 

299 

5 

44 

299 

9 

45 

299 

11 

46 

299 

13 

47 

299 

15 

48 

299 

17 

49 

299 

21 

50 

299 

3 0 

51 

299 

31 

52 

299 

35 

53 

299 

37 

54 

299 

42 

5 5 

299 

43 

5 6 

299 

47 

57 

299 

49 

53 

299 

53 

59 ■ 

299 

55 

cO 

299 

53 

61 

299 

62 


Benefit 
rate 


rate ‘-'Pper Lcv;er 

(Rs/ Hect~m) bound bcund 


- xcu 


Remark 



C 

0 

o 

o 

o 

o 


6 

6 

6 

Reservo i 

350 

350 

350 

initial 

2 

2 

2 

Conditio! 

60 

60 

60 


3 

3 

3 


3 

3 

3 


9 

9 

9 


a 

8 

3 


1 

1 

1 


5 

5 

5 


12 

12 

12 




n c, ' 


Table 

3.-4 

( 

\ m ^ 

ntd, ) 





1 

2 

1^- 

4 

5 

6 

n 

1 

8 

62 

293 

63 

0 

SO- 

30 

SO- 


63 

295 

64 

0 

SO 

3 0 

SO 


64 

299 

66 

f; 

3 

3 

3 


65 

295 

£5 

'■ 

z 

2 

2 


6 6 

259 

71 

C 

9 

o 

9 


67 

295 

76 

o 

R 

8 

8 

reservoir 

6 8 

25'' 

79 

c 

45 

45 

45 

condition 

65 

41 8 ’5^ 

a 1 

0 

20 

2C 

2 0 


7C 

299 

33 


a 

3 

P 


71 

2 59 

Q« 


10 

1 w 

10 


72 

295 

39 

c 

2 

2 

2 


n 

2 5 b' 

9C 

i../ 


OJ' 

3 


7 4 

295 

91 

's. - 

- n 

X 

1C 

10 


75 

299 

93 

T’ 

2 

2 

2 


7 6 

255 

3i; 

r-^ 

w 

2 

2 

2 


77 

1 

3 


200C 

G 

1183 


73 

2 


0 

lOQ 

C 

44 


79 

'k 

7 

c 

2100 

0 

1016 


80 

5 

7 

0 

300 

0 

37 3 


31 


S 

w 

25 CO 

0 

1394 


32 

' 

y 

c' 

G 

ICO 

0 

65 


33 

11 

15 

C 

ICO 

L- 

35 

Regul ated 

84 

13 

15 

0 

510 

C 

286 

f lows 

85 

15 

17 

c 

65 0 

0 

321 


35 

17 

IS 

0 

20C0 

c 

210 


87 

19 

8 

0 

1 6C0 

0 

0 


88 

3 

21 

0 

4400 

0 

1179 


39 

21 

24 

0 

4600 

0 

■2978 


9C 

'3 n 

33 

0 

180 

0 

56 


91 

31 

33 

0 

90 

c 

59 




159 


Table B-4 (...Contd.) 


1 

2 


4 

5 

6 

7 

92 

33 

40 

C 

320 

n 

115 

93 

35 

39 


1250 

c 

619 

9 4 

3 *7 

-J / 

39 

0 

5 0 

c 

30 

95 

39 

4C 

C 

1450 

c 

399 

96 

40 

41 

c 

isoc 


1C14 

97 

4di 

44 

o 

3CC 

r 

614 

93 

o ~ 

44 

C 

65 

r'-H 

49 

9S 

4 

41 

o 

67 C 

0 

543 

ICC 

41 

5 2 

c 

20'2C 

c 

1457 

1 01 

47 

51 

0 

33 

o 

19 

,102 

49 

5 1 

0 

135 

0 

70 

1C3 

51 

,J 45 

r-, 

lec 

c 

39 

104 

52 

24 

0 

25 CC 


15 46 

105 

2 7- 

61 

C' 

6300 

r 

1965 

106 

5 5 

57 

w 

120 

C 

12 

107 

55 

5 7 

c 

2C 

c 

4 

103 

5 7 

6C 

0 

19C 

n 

w 

16 

ICS 

5 3 

60 


19G 

0 

107 

110 

60 

61 

0 

205 

0 

123 

111 

61 

74 

o 

7250 

0 

2 CSS 

112 

62 

63 

0 

800 

c 

678 

113 

63 

67 

c 

320 

c 

769 

114 

66 

67 

0 

30 

c 

21 

115 

67 

59 

0 

7 40 

0 

640 

116 

69 

72 

0 

1230 

0 

1050 

117 

71 

72 

c 

35 

c 

11 

118 

72 

74 

0 

2200 

c ' 

1061 

119 

74 

78 

0 

7030 

0 

4509 

120 

76 

78 

0 

45 

0 

1 

121 

7 7 

88 

c 

9225 

0 

4510 

122 

79 

81 

0 

60 

0 

5 

123 

O 1 

ci X 

85 

0 

155 

0 

41 



Table T-4 (...Contd.) 


1 

2 

3 

4 

5 

6 

7 

8 

124 

33 

35 

T' 

25 

0 

6 


125 

G'5 

86 

o 

18C 

c 

47 


126 

86 

8S 

0 

530 

c 

195 


127 

: J O 

c c 


1 025 0 

0 

4895 


12 3 

po 

93 

U 

6C 

/-S 

27 

Regulated 

129 

P* 

98 

0 

SC 

0 

14 

flows 

130 

93 

97 

C 

65 

f- 

25 


131 

95 

99 

0 

25 

0 

6 


132 

97 

93 

0 

1 C'5 

r.; 

31 


133 

98 

99 

•s_>' 

210 

c 

45 


134 

99 

1 00 


1243C 

0 

5 090 


135 

1 

4 

45 

15 3 

c 

153 


136 

2 

4 

45 

17 

c 

17 


137 

5 

6 

45 

49 

n 

49 


133 

9 

10 

45 

16 

0 

16 


139 

11 

12 

45 

30 

0 

30 


140 

13 

14 

45 

IOC 

Q 

IOC 


141 

16 

14 

45 

0 

Q 

0 


142 

16 

13 

45 

35 

C 

35 

Irrigation 

143 

17 

13 

45 

22C 

220 

1 20 

f lows 

144 

17 

25 

43 

626 

0 

626 


] 45 

19 

6 

45 

210 

210 

21 C 


146 

g 

20 

45 

280 

280 

280 


147 

21 

22 

45 

650 

650 

650 


143 

21 

25 

45 

440 

0 

440 


149 

30 

34 

45 

70 

c 

70 


150 

3 1 

32 

45 

16 

16 

16 


151 

35 

34 

45 

500 

500 

500 


152 

35 

54 

45 

300 

0 

300 


153 

37 

38 

45 

6 

0 

6 






Contd. ) 


"able 3-4 (... 


1 

2 

3 

4 

5 

6 

7 

154 

44 

45 

45 

120 

120 

120 

155 

41 

46 

45 

100 

IOC 

ICG 

15 6 

A J 

43 

45 

Q 

0 

9 

1 7 

49 

50 

45 

20 

0 

20 

153 

24 

29 

45 

1000 

1000 

IGOO 

159 

- 

54 

45 

£C 

c 

60 

ISC 

55 


45 

3 

r"' 

3 

151 

53 

5 9 

45 

n 

f y 

12 

12 

162 

62 

25 

45 

3 4 

G 

34 

1 

55 

2 S 

53:7 

60 

60 

60 

16 4 

64 

65 

45 

2 C' 

c 

2C 

1 65 

64 

S2 


10 

c 

10 

165 

"'5 

65 

45 

7 

0 

7 

167 

67 

63 


150 

1 5 0 

150 

168 

69 

7 C 


20 

0 

7 r". 

169 

71 

7 75 
* 

45 

9 

2 C 

20 

170 

74 

75 

‘'ir 

550 

■K UZ 
*0 

65 ^ 

171 

7 6 

77 

/ 7 
s 

•3 

c 

3 

172 

81 

82 

45 

15 

15 

15 

173 

33 

84 

43 

15 

15 

. 15 

17 4 

86 

O’J 

45 

4C 

0 

40 

175 

pc 

92 

45 

5 

C 

5 

176 

90 

92 

45 

32 

0 

32 

177 

91 

92 

45 

20 

fy. 

20 

17 8 

91 

96 

45 

C 

0 

0 

179 

93 

94 

45 

16 

'w' 

16 

ISO 

95 

96 

45 

4 

0 

4 

181 

21 

28 

0 

1000 

1000 

1000 

132 

28 

74 

96 

1000 

1000 

1000 

183 

35 

3 6 

0 

250 

250 

250 




16 


Tatle B-4 (... Contd. ) 


1 

2 

3 

4 

5 

6 

7 

3 

134 

^ c> 

3 9 

96 

250 

250 

250 

Power supply 

1 85 

21 

23 

lOCO 

15 

15 

15 


-1 if- 

X w C 

21 

27 

ICCC 

15 

15 

15 

M & I aupply 

137 

35 

1:! 0 

icoo 

40 

w 

40 


123 

79 

SC 

lOCO 

/ 1 

* 

A 1 

■■'X 

43 

1 

1 39 

1 

151 

0 

145 

145 

145 


1 9 C 

2 

152 

o 

6 

6 

4 


191 

3 

153 

c 

1 33 

133 

93 


192 

Q 

159 

sr- 

19 

19 

9 


193 

11 

161 

0 

28 

28 

15 


194 

1 3 

163 

G 

27 

8 

19 


195 

16 

166 

c 

2 0 

20 

2C 


196 

17 

167 

'O 

9 6 

40C 

5 60 


197 

21 

171 


1140 

6CC 

540 


198 

3C 

130 

0 

74 

74 

34 


199 

31 

181 

0 

11 

6 

5 


2 00 

35 

1 QC 

X O— 

G 

5 41 

350 

191 

Reservoir 

201 

202 

37 

4-2 

137 

192 

p. 

C 

C 

kj 

122 

2 

60 

4 

62 

inif ,al 
condition 

203 

43 

193 

0 

D 

3 

5 


204 

47 

197 

c 

7 

7 

4 


9 

49 

199 

0 

29 

29 

20 


206 

53 

203 

0 

31 

31 

23 


207 

5 5 

2 ud 

0 

4 

4 

3 


2C3 

53 

203 

0 

10 

5 

D 


209 

G2 

212 

0 

40 

12 

2S 


210 

63 

213 

c 

64 

30 

34 


211 

64 

214 

ns 

K.. . 

60 

30 

30 


212 

66 

216 

n 

7 

3 

4 


213 

69 

219 

0 

2 

■ 2 


I 

214 

71 

221 

0 

13 

9 

9 








163 


“bis 3-4 Ccntd. } 


1 

2 

3 

4 

5 

6 

7 

s 

215 

1 

22 6 

w 

24 

S 



216 

79 

2 2 '1 

c 

82 

45 

37 


217 

■31 

•7 7 1 
.j- 

c 

7 c 

2C 

19 


21 - 

83 

2:3 



S 

9 


219 

35 

23 6 


42 

42 

32 

Reservoir 

22C 

39 

23 9 

'C' 

5 

5 

3 

initial 

22i 

9 C: 

24C 

n 

22 

22 

14 

conditions 

22 2 

91 

241 

c 

35 

20 

10 


223 

G 1 

2 43 

5' 

S 

s 

6 


224 

93 

245 

G 

C 

0 

0 


m 

• 

m 

m 


• 

* 

• 

Inflows 

m- 

•■ 

m 

♦ 

* 


. 

(Table 5,12) 

265 

151 

153 

C 

132 

132 

40 


266 

152 

153 

c 

7 

7 

0 


267 

153 

157 

c 

139 

139 

40 


268 

155 

157 

0 

55 

55 

90 


269 

157 

158 

c 

194 

194 

130 


27 C 

159 

153 

0 

20 

20 

3 


271 

161 

165 

0 

ICO 

G 

5 


27 2 

163 

165 

c 

510 

C 

0 


273 

165 

167 

0 

65 0 

0 

c. 


274 

167 

169 

0 

2000 

G 

125 


275 

169 

157 

c 

1 600 


cC 

Regnl ated 
flows 

276 

15 3 

7 r* -1 

X # X 

0 

4400 

0 

jO, 

u 


277 

171 

174 

0 

4600 

O' 



273 

13C 

183 

G 

180 

0 

24 


279 

1 31 

133 

c 

90 

r- 

6 


2 8C 

1S3 

190 

c 

320 

. 0 , ' 

30 


281 

1 85 

139 

c 

1250 

0 

0 


282 

187 

139 

0 

150 

0 

9 




Contd, ) 


Table 3—4 


123 45 67 8 


2 33 

1C9 

190 

0 

1450 

234 

190 

191 

0 

1800 

285 

192 

194 

0 

500 

236 

193 

154 

c 

65 

2 37 

194 

191 

V.” 

670 

238 

191 

2 02 

c 

2020 

2 39 

197 

201 

C' 

38 

29 C 

199 

201 

c 

135 

29i 

201 

202 

0 

ISC 

292 

202 

174 

n 

2500 

293 

17 4 

211 

0 

6300 

294 

203 

207 

c 

120 

2S5 

205 

0 

0 

20 

29 5 

207 

2 1 0 

0 

19G 

297 

2 '-'d 

21 C 

0 

190 

29S 

210 

211 

0 

205 

299 

211 

224 

/> 

W' 

7250 

3 OC 

J, -C-- 

213 

c 

300 

3 01 

213 

.b b. « 

0 

S20 

3C2 

* 'C^ 

217 

c 

3 0 

3 03 

217 

219 

r\ 

c 

74C 

3 04 

219 

222 

0 

1230 

305 

221 

222 

c 


3C6 

222 

224 

r'- 

2 2 CO 

307 

224 

22S 

c 

9030 

3 03 

226 

223 

0 

45 

309 

? 2 9 

23 3 

n 

9225 

310 

229 

231 

c 

60 

311 

231 

2 ■’■5 

c 

155 

312 

233 

235 

0 

25 


C 

259 


0 

239 


0 

106 


0 

11 


n. 

*>w<i 

45 


0 

23 4 


c 

1 


0 

21 


c , 

22 


0 

256 


0 

0 


c 

3 


c 

3 

Regulated 

0 

11 

flOV7S 

0 

3 


0 

14 


0 

14 


0 

p p 


c 

109 


0 

A 


r, 

'w 

0 


o 

c 



0 



l' 


u 

504 


G 

14 


0 

513 


0 

5 


c 

19 


7 

7 

V.’ 




16 


Table 3-4 (... Contd. ) 


± 

2 

3 

4 

5 

6 

7 

5 

313 

235 

236 

0 

IBO 

0 

19 


314 

23 5 

23 3 

0 

530 

C 

41 


313 

2 3 8 

249 

.0 

10250 

0 

579 


316 

239 

2 43 

0 

60 

c 

4 


317 

24C 

243 

0 

60 

c 

3 

Regtil ated 
flows 

313 

243 

247 

0 

65 

0 

4 


319 

245 

247 

0 

25 

c 

3 


320 

247 

243 

u 

105 

pi 

7 


3 21 

O •* D 

4-5 

Z 49 

r. 

w 

21C 

n 

IG 


322 

249 

250 

U 

1 2 43 0 

6C'C 

6 z 


3 23 

151 

15 4 


105 

C 

105 


3 24 

15 2 

J. b ‘-x- 

A 

1C 

c 

10 


3 25 

155 

156 

45 

33 

c 

3 3 


326 

15 9 

16C 


16 

0 

16 


327 

161 

162 

45 

20 

0 

20 


323 

153 

164 

45 

44 

c 

4 A 

Ttb 


329 

166 

1 64 

/I r; 

‘•jT'wI 

20 

n 



33C 

166 

16.3 

A, ^ 

20 

C' 

35 


331 

1 67 

163 

45 

260 

140 

140 


332 

167 

175 

45 

330 


3 40 


333 

1 69 

15 6 

45 

210 

90 

9 0 

Irrigation 

33 4 

153 

170 

45 

280 

220 

22C 

flows 

335 

171 

172 

45 

65 C 

200 

519 


336 

171 

176 

45 

22S 

0 

0 


•-S n r-7 
/ 

ISO 

134 

45 

30 

0 

30 


333 

1 01 

1S2 

45 

16 

9 

9 


33S 

135 

134 

45 

5 CO 

50 

101 


3 4 0 

1S5 

204 

45 

200 

0 

■ Q 


341 

137 

138 

45 

4 

■ 0 

4 


3 4 2 

194 

195 

45 

120 

72 

12 




Table B-4 ( 


Cental. ) 


1 

2 

3 

4 

5 

6 

7 

8 ; 

3 43 

191 

196 

45 

100 

100 

100 


3 44 

197 

198 

45 

6 

0 

6 


3 45 

199 

200 

45 

9 


9 


3 46 

17 4 

179 

45 

1000 

2-3 

406 


3 47 

203 

204 

45 

35 

c 

35 


3 43 

2 '05 

2C6 

45 


n-. 

2 


3 49 

208 

2C9 


12 

p 

S 


35 0 

2 1, 2 

i / o 

45 

23 

c 

23 


351 

213 

176 

45 

A n 

40 

40 


352 

2 1 4 

215 

45 

50 

0 

20 


353 

214 

23 2 

if K 

30 

c 

10 

Irrigation 

354 

216 

215 

45 

3 

3 

5 

f lov/s 

355 

217 

21G 

4 o 

150 

ICO 

110 


35 6 

215 

220 

45 

X 

c 

12 


3 57 

221 

223 

45 

20 

14 

14 


353 

224 

225 

45 

65 C 

Ti-OO 

106 


359 

225 

227 

45 

A 

P 

4 


3 60 

231 

232 

45 

25 

0 

5 


3 61 

233 

23 4 

45 

15 

12 

14 


362 

236 

237 

45 

20 

'■sj 

20 


3 62 

239 

242 

45 

c 

0 

c 


364 

240 

242 

45 

21 

0 

21 


3 65 

241 

242 

45 

20 

0 

20 


3 66 

241 

246 

45 

20 

p 

O 


367 

243 

244 

45 

9 

c 

9 


3 58^ 

245 

246 

45 

2 

G 

2 


359 

171 

17 8 

C 

lOOO 

1000 

1000 


37C 

1 7 

224 

96 

lOOO 

lOOO 

1000 

Pov;er 

371 

185 

186 

0 

250 

250 

25 0 

flOVJS 

372 

135 

189 

96 

250 

250 

250 




1 6 


ral:la 3~4 Contd, ) 


*? 

2 

3 

4 


c: 

7 

8 

573 

171 

173 

lOCO 

15 

1 " 

1 5 


37 4 

171 

17? 

ICCO 

lb 

15 

1 c, 

i .J 


375 

185 

230 

iOOC 

40 

40 

4C 

M 6 I flows 

376 

229 

230 

icoo 

42 

42 

42 


3/7 

151 

3 00 ^ 

w 

1 

ICO 

ICO 


373 

152 

3 CO 

W‘ 

2 

2 

2 


379 

15 5 

3 OG 

C 

9C 

90 

90 


3 SC 

159 

3 CO 

0 

10 

1C 

10 


3 31 

131 

3 CO 

w 

12 

12 

12 


3 32 

163 

2 00 

c 

o 

8 

S 


3 33 

ic5 

5 00 

c 

20 

20 

20 


3 34 

1 

3CC 


300 

3 CO 

300 


3 35 

171 

3 CO 

sT- 

6 CC 

6 -00 

600 


3S6 

1 SO 

3CC 

r'. 

O' 

4C 

40 

40 


3 37 


3 00 



6 

5 


3 83 

1 o tz 

3 CC 

n 

35 0 

350 

350 

Reservoir final 

3S9 

137 

O, 

n 

2 

2 


390 




2 

conditions 

192 

3 CO 

0 

60 

60 

60 


3S1 

1 ■*3 

3 00 

0 

-.J 

3 

3 


392 

157 

3 CO 

c 

3 

3 

3 


393 

199 

3 OG 

c 

9 


9 


394 

2 03 

3 CO 

G 

P 

8 

8 


395 

205 

3 00 

V--' 

1 

1 

1 


39 5 

20^' 

oo 

n 

- 

X 

c; 


397 

212 

3 00 

C 

12 

12 

12 


393 

213 

300 

L,' 

30 

30 

30 


33? 

214 

300 

0 

30 

30 

30 


40G 

216 

300 

0 

3 

3 

3 


401 

219 

3 00 

0 

2 

2 

2 


402 

221 

30C 

0 

9 

9 

9 



IS 


i able 

3-4 

C « • • 

Ccn 

td. ) 




1 

2 

3 

4 

5 

6 

7 

3 : 

4C3 

226 

SCO 

0 

s 

8 

8 


4G4 

229 

30C 

0 

45 

45 

45 


4 C5 

231 

1 i-.r^ 

ji 

0 

20 

20 

20 


4C6 

23 3 

300 

c 

3 

S 

8 

Final reser- 

<4: w t 

23 6 

3 00 

r. 

10 

10 

iO 

voir conditivcr 
(Rabi seascr.' 

■^r ‘~''C3 

239 

300 

r- 

2 

2 

2 


,n r**' 

24C 

00 

C 

3 

3 

3 


410 

241 

3 00 

C 

10 

10 

10 


411 

243 

1 nn 

0 

2 

2 

2 


412 

245 

*5 00 

C 

4 

4 

4 


413 

4 


C 

115 

0 

112 


• 

* 


• 


• 

0 


« 

• 



* 

m 

0 


* 




m 

m 

m 

Exit flows 


• 


* 

• 

m 

* 


* 

• 


m 


* 

# 


- 

• 


• 

• 

- 

m 


4S7 

J b. 1.. 

299 

C 

2073C 

0 

15 404 





169 





Table B , 

Value of 

5 



Arcs 

1-5 

6 

7^9 

10 

11 

12 


C 

45 

c 

45 

c 

45 

Arcs 

13 

14 

i- b 

± o 

17 

13 


0 

4 b 

f'-, 

45 

C 

45 

Arcs 

19 


21 

22 

23-2 4 

25 


C 



45 

0 

45 

Arcs 

2 6-27 


29-35 

3 6 

37 

3 2 


G 

45 

0 

45 

0 

45 

Arcs 

39-49 

5 0 

51-69 

IQ 

71-36 

37 

'S'- 

0 

45 

‘w’ 

45 

■ r 

45 

A.rcs 

88-9 C 

91-92 

93 

94 

95 

96 

cx- 

0 

45 

C 

45 

C 

45 

Arcs 

97-1 GO 

101-150 

151^152 

153 

154^155 

156-157 


C 

45 

C 

45 

. ■’ 

45 

Arcs 

153-159 

1 ow 

161 

162 

163 

164 

CV* 

0 

45 

Q 

45 

C 

45 

Arcs 

165-177 

172 

179--180 

131 

1 32-1 S4 

1 OS -136 


0 

45 

c ■ 

45 

0 

45 

Arcs 

137 

138 

139-197 

197 

199 

200 

o- 

0 

45 

w 

45 

f\ 

45 

Arcs 

201-203 

204 

2C6 

2C5 

207^219 

220221 


0 

45 

0 

45 


45 

Arcs 

222 

223 

224-226 

227 

2 2 ' : ■ 

229 


C 

45 

0 

45 


45 

Arcs 

230-232 

23 3 

234-236 

237 

238-240 

241-242 

O'' 

0 

45 

0 

45 

0 

45 

Arcs 

2 43 

244 

245 

246 

247-250 

251-293 

o- 

Arcs 

299-487 

45 

0 

, 45 


45 


<y~ O 




171 


1 

2 

7 

.4 

■~x 

5 

94 

-7 

39 

G 

2 7 , 

73 

IQ 

G r 

fG 

i 5 9 

9 c 

3C 

.■*, 1 


“ GG- 

1 / / 

42 

GG 

: 

orsn 

9 3 

43 

G . -I 

C- 

7 7 

o o 

44 

.1 

0 

67 0 

1 3'“ 3"'"' 

41 

3 


202 c 

"5, .r-: "I 

1 

47 

7 " 

c 

7 G, 

1C2 

Q 

b 

G, 

135 

1 C'3 

51 

G G 

G 

1 30 

104 

52 

2 4 

C 

^ jT^ r^.f^ 

1 03 

24 

G1 

r~’>, 

53GO 

1C5 

53 

7 


120 

1 G7 

55- 

53 

c 

20 

1C3 

57 

< G> 

kJ 

190 

109 

5 7 

6g' 

0 

190 

110 

60 

o X 

c 

205' 

111 

Ox 

7 4 

0 

G 7 G' 

112 

62 

53 

c 

3 CO 

1 1 3 

63 

67 

r) 

320 

114 


r-j 

C / 

c 

7 0 

115 

67 

5e 

n 

/ 4^-2^; 

i 1 S 

u :: 

7 7 

0 

1230 

117 

71 

* 

C 

35 

113 


e gX 

C 

2200 

119 

7 4 

7 3 

0 

903 0 

120 

7 6 

73 

0 

45 

121 

7 3 

O '' 7 . 

c 

9 2 25 

122 

y G 

■ 81 

0 

60 

123 

31 

. 85 

. c 

1C':; 
j. J -J 


7 


a 



-- 

0 


0 

G C G 

7 : 

4 1 0 


35 

W 

G G 7 

0 

773 


10 

C: 

Gf"-’ 

0 

50 


C 333 

n 

Regulated 
C 0 flows 

C 2 

0 2 



54 

0 

■ 60 

u 

734 


/ Gj 0 


~5r 

i. 

G J 4 

c 

■ 15 

c 

G 0 Cl 

c 

777 

C.' 

1 

C: 

7 0 4" 

4*'' 

O-.' 

7^3 O' 

r\ 




0 

2034 

C 

■ O' 

c 

26 



25 




17 2 


1 

z 

3 

o- 

5 

6 

7 :-| 

'53 

7 7 

3 2 

'■5 

2 

C; 

6 

15 4 

4* 

■ CT' 


120 

12C 

120 

1 5 5 

4-i 

• 1- b 

ob 

1 C2 

ICC 

ICC) 

1 K 



. !*• 




Jl, V 



VJ 

4 


9 ' 

1^7 


- ■ 

. - 








2',.. 

■ 0 

2C 

153 

'^5' ,/T 

■Z* *x 

47 

J 

1 ocr: 

iccc 

ICOC 

159 

53 

5 

’4x 

sc 


SO 

16C 

55 

rr )C.. 



0 

3 

161 

7 p 

5 9 

/'7. 

12 

12 

12 

1 62 

62 

25 

/ ” 

34 

r- 

3 4 

163 

63 

2 5 

4j 

60 

so 

5C 

154 

64 

65 

‘x3 

35 

r*-. 

20 

155 

64 

32 


35 

c 

1C 

156 

< c 

DO 

25 

■1 n 

O 

7 

5 

157 

67 

6S 

45 

150 

150 

150 

165 

69 

7 p, 

45 

20 

C 

20 

1 -c.g 

71 

^ 7 

/ .J 

45 

20 

2C 

20 

170 

74 

75 

45 

65 C 

55 C 

55 C 

171 

76 

77 

45 

S 


n . 

172 

31 

O 9 

40 

25 

or, 

1 5 

173 

77 

p 

45 

15 

15 

15 

17 4 

36 



40 


40 

175 

39 

7 7 

45 


p. 

5 

17 5 

90 

92 

45 

32 

C 

3 2 

177 

O 1 

9 2 

45 

2C 

0 

20 : 

17 1 

91 

96 

45 

20 

b' 

o 

J- / 

93 

O 

45 

1 7 

0’ 

16 

“t O'"' 

95 

95 

45 

4 

0 

4 ;/ 

131 

21 

23 

C 

icco 

■IC-OC 

l-'OO 

132 

23 

7^^^ . 

96 

. 1000 

■lOCC.- 

1000 



173 


Table B-6 ( . . . Contd. ) 


153 

37 

38 

45 

6 

0 

6 

15 4 

44 

45 

45 

120 

120 

120 

155 

41 

46 

45 

100 

100 

100 

156 

47 ■ 


45 


0 

9 

157 

49 

5 0 

45 

20 

0 

20 

158 

24 

o o 

45 

1000 

1000 

1000 

159 

53 

54 

45 

60 

0 

SO 

160 

55 

55 

45 

3 

0 

3 

161 

5 8 

59 

45 

12 

12 

12 

162 

62 

25 

45 

34 

0 

3 4 

163 

63 

26 

45 

60 

60 

60 

164 

64 

65 

45 

35 

0 , 

20 

165 

■ 64 

32 

45 

35 

c 

10 

166 

66 

65 

45 

3 

8 

5 

167 

67 

68 

45 

150 

150 

150 

168 

69 

70 

45 

20 

C 

20 

169 

71 

73 

45 

20 

20 

20 

170 

74 

75 

45 

650 

65 0 

65 0 

171 

76 

77 

45 

8 

0 

8 

172 

81 

82 

45 

25 

25 

15 

173 

83 

84 

45 

15 

15 

15 

174 

S6 

37 

45 

40 

r-. 

40 

175 

89 

92 

45 

O 

0 

5 

176 

90 

92 

45 

32 

0 

32 

177 

91 

9 2 

45 

20 

0 

20 

17 3 

91 

96 

45 

20 

0 

0 

179 

93 

9 4 

45 

IS 

0 

16 

1 80 

95 

96 

45 

4 

0 

4 

181 

21 

28 

0 

1 OCO 

1 000 

lOOO 

132 

28 

74 

96 

1000 

1000 

1000 



TaiDle E-S 


( . . . Contd. ) 


1 

2 

3 

4 

5 

o 

/ 

8 

183 

35 

36 

0 

250 

25 0 

25 0 


184 

36 

3S 

96 

250 

25 0 

250 


135 

21 

23 

1 OCO 

15 

15 

15 


136 

21 

27 

lOCO 

15 

15 

15 


1 37 

35 

30 

1000 

40 

40 

40 


1 3 3 

79 

80 

1000 

43 

43 

43 


139 

1 

J_ 

151 

0 

145 

145 

100 

Reservoir 

. 

♦ 

* 

• 

* 

* 


initial 

. 

• 

* 


• 


* 

condition 

2 2 4 

95 

2 45 

0 

4 

4 

4 

rabi saaso: 

225 

299 

151 

0 

110 

110 

75 

Inf lox'.'s 

<■ 


. 

o 

• 


. 

raloi 

• 

• 

- 


* 

- 

• 

season 

264 

299 

2 49 

0 

37 

37 

21 


265 

151 

153 

r-, 

V./' 

2000 

2 OCO 

15 


266 

152 

153 

r-. 

100 

ICC 

0 


267 

153 

157 

0 

2100 

2100"^ 

15 


268 

155 

157 

0 

8 00 

80C 

32 


269 

157 

158 

c 

25 00 

25 00 

97 


27 C 

159 

15 3 

0 

100 

100 

G 


271 

161 

1 CK 

j, 

o 

IOC 

0 

2 


272 

163 

165 

0 

510 

c 

0 


273 

165 

167 

0 

65 0 

0 

2 


274 

167 

169 

c 

2 OCO 

0 

290 


275 

169 

15 3 

0 

1600 

c 

200 


276 

153 

171 

c 

4400 

c 

20G 


277 

171 

174 

C) 

4600 

o 

0 


278 

180 

133 

c 

1 GO 

0 

13 


279 

131 

133 

0 

90 

G 

0 


280 

133 

190 

0 

320 

o 

13 


281 

185 

1S9 

c 

1250 

6 

0 




175 


Z'd ^ 
2S9 

290 

291 
29 2 

293 

294 


191 

197 

199 

201 

202 

174 

203 


202 

201 

201 

202 

174 

2 'I 1 


C 

0 

o 

0 


Table 

iD-6 1 

L . » Conta 

• / 



1 

2 

3 

4 

5 

6 

232 

137 

139 

0 

150 

r^. 

O' 

233 

139 

1 Q O 







f 5 

145G 

c 

2S4 

190 

191 

r-. 

1300 

0 

235 

192 

194 

C 

SCO 

C 

236 

193' 

194 

Ly 

55 

0 

z 

194 

191 

C 

670 

0 


2020 
18 
135 
ISC 
25 00 
6000 
120 


O 

O 

O 

o 

c 

o 

0 


6 

25 6 
274 
95 
10 
33 
2 07 
0 

18 
18 
225 
O 
n 


295 

205 

.207 








■o 

20 

0 

1 

296 

207 

9 1 r* 







O -1. w' 


190 

L. 

1 

297 

208 

O 1 ri 

-L. V.7 

c 

190 

0 

'0 

293 

210 

0 17 







•eO _L J_ 

C' 

205 


1 

299 

211 

22 4 

0 

7250 

0 

i. 

3 00 

212 

217 

0 

300 

0 

65 

301 

213 

217 

c 

320 

0 

SO 

302 

21 S 

217 

0 

30 

c 

2 

3 03 

217 

219 

0 

7 40 

0 

0 

3 04 

219 

2 22 

o 

123 0 

r> 


3 05 

221 

222 

'Z 

7'- 


0 

3 05 

222 

224 


22 00 

0 

0 

3 07 

224 

223 

c 

9030 

c 

480 

308 

226 

22 3 

c 

45 

o 

12 

3 09 

223 

23 3 

0 

9225 

0 

492 

310 

229 

231 

V-' 

60 

0 

1 

311 

231 

235 

0 

155 

0 

.10 

312 

233 

235 

0 

25 

c 

n 



Table B-6 ( 


Contd. ) 





1 

2 

3 

313 

235 

235 

314 

236 

233 

315 

23 3 

249 

315 

239 

243 

317 

240 

243 

318 

243 

247 

319 

245 

247 

32 0 

247 

243 

3 21 

243 

249 

322 

249 

250 

3 23 

151 

15 4 

3 24 

152 

154 

325 

155 

1 45 

3 2 6 

159 

15C 

327 

ISl 

162 

3 2S 

163 

1 

329 

166 

1 ^^4 

330 

1 66 

i 6S 

331 

167 

163 

33 2 

167 

175 


/■ C 

45 

45 


20 
20 
2 SC 
3P''' 


O 

O 

140 


333 

169 

156 

45 

210 

90 

33 4 

153 

17C' 

45 

23C 

220 

335 

171 

172 

45 

65 0 

535 

33 6 

171 

1 *7 . 0, ■ 

45 

223 


337 

180 

184 

45 

30 

0 

333 

181 

109 






16 

9 

339 

. 135 



5 00 



-L 

4o 

4CC 

3 40 

185 

2 C4 

45 

200 

C 

3 41 

137 

133 

45 



342 

194 

195 

45 

120 

72 


4 

5 

6 

7 

O 

1 30 

0 

10 

0 

530 

0 

29 

0 

1C25C 

c 

534 

0 

60 

c 

1 

G 

60 

c 





u 

0 

65 

c 

0 

0 

25 


1 

0 

105 

0 

1 

0 

210 

0 

1 

c 

12430 

0 

55 6 

45 

105 

c 

105 

45 

1C 

G 

3 

45 

33 

c 

3 3 

45 

IS 

c 

16 

45 

20 

0 

CC 

45 

44 

0 

40-1 


0 

25 

115 

397 

90 

97 

20 

n 


72 


177 


Table B-6 (...Contd.) 


1 

2 

3 

^■3 

5 

6 

7 

3 43 

191 

196 

<45 

100 

100 

ICC 

3 44 

197 

1 O O 

1 b' o 

45 

6 

0 

6 

3 45 

199 

200 

45 

9 

c 

9 

346 

174 

179 

‘ tD 

1 COO 

1000 

335 

3 47 

203 

2C4 

45 

35 

0 

33 

3 43 

205 

206 

45 

2 

0 

2 

349 

20,3 

205 

45 

12 

s 

7 

350 

212 

175 

45 

23 

c 

23 

351 

213 

176 

45 

50 

40 

40 

352 

214 

215 

4-5 

3 0 

c 

20 

353 

214 

232 

A" 

30 


C 

35 4 

216 

215 

45 

0 

0 

6 

r 

355 

217 

21? 

-'.5 

150 

•82 

82 

356 

219 

220 

45 

16 

0 

5 

357 

221 

223 

45 

20 

14 

12 

353 

22 4 

225 

45 

65 C 

550 

550 

359 

226 

227 

45 

/! 

0 

,4 

‘‘X 

3 50 

231 

232 

45 

25 

15 

5 

361 

233 

23 4 

45 

15 

15 

14 

362 

236 

237 

45 

20 

0 

20 

3 63 

239 

242 

45 

4 

0 

5 

354 

240 

242 

45 

21 

c 

20 

365 

241 

242 

45 

30 

0 

20 

3 66 

241 

246 

45 

30 

0 

0 

367 

243 

244 

45 

9 

c 

8 

368 

245 

246 

45 

2 

c ■ 

.2 

369 

171 

178 

0 

1000 

1 000 

ICCO 

370 

178 

224 

96 

1000 

ICCC ' 

IGCC 

371 

135 

185 

0 

250 

25 0 

250 

372 

186 

189 

96 

250 

25 0 

250 



17 S 


Table 

3-6 ( 

... Gontd. 

) 




* 

1 

2 

3 

4 

5 

6 

1 

3 

373 

111 

173 

1000 

'"'15 

15 

15 


374 

111 

177 

lOCO 

15 

15 

15 


375 

X 

23 0 

1000 

40 

40 

40 


376 

229 

23 0 

1000 

4: 

42 

42 


377 

151 

30C 

o 

1 ''.'O 

100 

♦ 

100 

Exit 

flows 

4B7 

3 00 

299 

c 

20731 

• 

c 

16704 




179 

Table B-7 
Value of i?v- 


Arcs 

1-3 

4 

5 

6 

7-9 

lO 


0 

45 

C 

45 

0 

45 

Arcs 

11 

12 

13 

14 

15 

16 


0 

,n C, 

C 

45 

c 

45 

Arcs 

17 

11 

1 3 

20 

21 

2 2 


0 

45 

0 

45 

0 

45 

Arcs 

23-24 

25 

2S~27 

28 

29-35 

36 


0 

45 

w 

45 

G 

45 

Arcs 

37 

■3 Q 

.3 C 

39-47 

43 

49 

50 

<3^ 

0 

45 

0 

45 

C 

45 

Arcs 

51-53 

5 4- 

55 

5 6—69 

70 

71-76 

c?'. 

0 

45 

0 

45 

0 

45 

Arcs 

77 

7S-S6 

37 

88-90 

91-92 

93 

O-^ 

0 

45 

0 

45 

Q 

45 

Arcs 

94 

95 

96 

97-100 

101-150 

151-152 


45 

' 0 

45 

0 

45 

0 

Arcs 

153 

154-155 

156-157 

158-159 

160 

161 


45 

0 

45 

0 

45 

0 

Arcs 

162 

163 

164 

165-177 

178 

179-130 


45 

0 

45 

0 

45 

0 

Arcs 

181 

132-134 

185-136 

137 

138 

139-197 


45 

0 

45 

0 

45 

0 

Arcs 

198 

199 

200 

201-203 

204 

205 


45 

0 

45 

0 

45 

0 

Arcs 

206 

207-219 

220-221 

222 

223 

224-226 


45 

Kj 

45 

0 

45 

0 

Arcs 

227 

223 

229 

230-232 

233 

234-236 


45 

0 

45 

0 

45 

0 

Arcs 

237 

233-240 

241-242 

243 

244 

245 


45 

0 

45 

0 

45 

0 

Arcs 

246 

247-250 

251-298 

299-487 



12 -^' 

45 

0 

'45 

0 





RSPEREIICES 


Becker, L. and Yeh (1974) . Optimal timing, sequencing a 
sizing of multiple reservoir suxface v/ater supply 
facilities. Hater Resour, Res. 10(l), 57-65. 

Butcher, W. S., Hall, 71. A., and Kaimes, Y.Y. (1969). Dynai 
prograrroning for the optimal sequencing of water suppi 
projects. 7Jater Rosour. Res. 5, 1196-1207. 

Buras, I-I. (1972). Scientific allocation of water reso’crc! 
Amer, Elsevier, Mow Yorlc. 

Burton, R.M. (1969) . A model for project selection in 

developing countries, 36th national ORSA meeting, Mi-: 
Florida. 

B7JSS3 (1974). Bangalore Hater Supply and Sewerage Board, 
Bangalore Hater Supply Reg’^ilations, Karnataka Gazett*; 
(extraordinary), March 12, 1974. 

Charnes, A., 7J.W- Cooper, and G. H. Syraonds (1958). Cost 
horizon and certainty equivalents on approach to 
stochastic programming of heating oil. Management 
Science, 4(3), 235-263. 

Chow, V.T. , and kJindsor, S. (1972). Multi reseroir optim 
tion model. Proc. A. S.C.E., 98, Hyd 10, 1827-1842. ; 

CPF (1972) . Report of the Cair/ery fact finding Comrrdttd 
Ministry of Irrigation & Pov/er, Government of India, i 

Chun, R.Y.D. , et al. (1964). Ground li'ater management for 
nations f uture-optimi.im conj motive operation of groun 
water basins. Proc. A. S.C.E., 90, Hyd. 4. 

CW&PC (1955). A short note on increasing government incoj 
jDy raising charges for water. Central Water and Power 
Commission, Ministry of Irrigation & Power, New Delhi 

Dracup, J.A. (1966). Optimum use of ground-water and | 
surface water system. A parametric linear programming 
approach. Tech. Rep. 6-24, Hyd. Lab., Univ. Calif., [ 
Berkeley, Calif. 


180 



isi 


Ecstein, O. (1965). Water resources development. Harvard Univ, 
Press, Cambridge, Mass., 1965. 

Eulkixrson, D.R. (1961). An out of Kilter method for minimal 
cost flows, SIAM, J. A.ppl. Hath., 18-27. 

Hall, W.A. and Dracup, J.A. (1970'). Water resources system 
engineering, McGraw-Hill, New York. 

Hadley, G. (1964) . Eonlinear and dynamic programming, Addison- 
V-Jesley, Reading, Mass., 1964. 

Harbeck, G.E., Jr. (1966). Annual variability of lake 

evaporation, Int. Ass. Sci. Hydrol. (Garda), Pu!d1.70(i), 
294-303. 

Hatfield & Pattison (1970) . The conjunctive use of water 

resources. Proc. of Ground Water Symp., The Univ. of New 
South Wales, Australia. Ed. K. K. Watson. 

Himmelblau, D.M. (1973). The optimal expansion of a water 

resources system. Final Technical Report, Dept, of Chsn. 
Eng. , Univ. Texas, Austin, Texas. 

Hirshliefer et al. (1960). Water supply; Economics Technology 
and policy. Univ. Chicago Press, Chicago. 

Hufschmidt, M.M. and Fiering, M.B. (1966). Simulation technique 
for design of water resource systems. Harvard Univ. Press, 
1966. 

KEB (1974). Karnataka Electricity Board. Electric Power Tariff 

Kimi, N.J. (1968). Linear programming with random requirements i 
Utah Water Res. Lab. Rep. PRWG 42- 2T, Utah State Univ., Loc 

Koening, L. (1963). Economics of ground water utilization. 

Amer. Water Works Assoc. Jour. 55 (l), 59-66. 

Loucks (1969) . Stochastic model for analyzing river basins. i 
Water Resour, and Marine Sci. Center, Cornell Univ., 

Ithaca, New York. 

Loucks (1968). Some comments on the optimization methods for ! 
branching multistage water resources systems. Water Resourj 
Res. 4(2), 447-450. 

Maass, et al. (1962). Design of water resources .systems. 
Harvard Univ. Press, Cambridge, Mass. 



182 


Manne, A.S., (1962), Product mix alternatives. International 

Eco . Review, 3 ( 1 ’) . 

Manne, A.S., (1967). Investments for capacity expansions, size, 

location and time phasing, M.I.T. press, Carcibridge, Mass. : 

Mao, C.T. (1969). Quantitative analysis of financial decision! 
MacMillan, London. 

Marglin, S.A. (1963). Approaches to dynamic investaent plannii 
contribution to economic analysis, Xl-CIX. ‘'msterdam, North- 
Holland PuId. Co. 

Milligan, H. (1970) . Optimizing conjunctive use of ground 

water and surface v/ater. PRVJG 42- 4T, Utah Water Research 
Laboratory. Utah State University, Logan, Utah. 

Morin, T.L. (1973). Optimal capacity expansion projects, 

Proc. ASCE 99. Hyd 9, 1605-1612. 

Morin, T.L. and Esogbue, M.O. (1971). Some efficient dynamic 
programming algorithms for optimal sequencing and schedulij 
of water supply projects. Water Res. Res. 7(3). 

Morin, T.L. (1973). Pathologg^ of a dynamic programming 

sequencing algorithm. Water Resour. Res. 9(5) 1178-35. 

Mullvihill, M.H. and Dracup, J.A. (1974). Optimal timing and : 
sizing of a conjunctive urban supply and waste water syster 
with non-linear programming. Water Resour. Res. 10(2), 170| 

Mayak, S.C. and Arora, S.R. (1971). Optiraal capacities of 
multireservoir system using the linear decision rule. 

Water Resour. Research. 7(3), 485-498. 

Nayak, S.C. and Arora, S.R. (1970). A separable prograi'nming 
approach to’ the determination of optimal capacities. 

Water Resour. Res. 9(5), 1166-1135. 

NCxAER (1959). Criteria for fixation of water rates. National 
Council of Applied Economic Research, Asia Publishing 
House, Nevr Delhi. 

O'Lagaire, D.T. and Himmelbleau, D.M. (1971). Optimal 

capital investment in the expansion of an existing water 
resources system. Water Resour. Bull. 7(6), 1194-1209, 

Orlob, G.T, (1970). Planning the Texas water system, Annu. 
Water P.esour. Res. Conf. 5th, Washington, D.C., 197 0, 
pp 171-136. 



1S3 


Raghava Rao, et al. (1969). An estimate of the ground water 
potential of India - First approximations. Soil & Water 
Management Symp.^ Hissar, March 1969. 

Revelle, C. et al. (1969). The linear decision rule in reserve 
management and design, 1 , Development of Stochastic model. 
Water Resour. Res. 5(4). 767-777. 

Rogers, P. and D. V. Smith (1970). The integrated use of grourt 
and surface waters in irrigation project planning. 

Amer. Jo’-ir. Agrl. Eco. 52 (i), 13-24. ; 

TDinJB (1968) . The Texas water plan. Texas Water Development 
Board, Austin. 

Tocher (1963). The art of simulation. The English University 
Press, London, 1963. 

Trott, W. and Yeh, W.I'J.G. (1973). Optimization of multiple 
reservoir system. Proc. ASCS Hyd 10. 1865. 

T^^son, H.U. and Weber, E.M. (1964) . Cbmputar simulation of 
ground water basins, Proc. ASCE 90, Hy 4, 

U. P. Irrigation Commission (1959). Fixation of water rates - 
Commission report. 

Veissman (1974). Slkhorn river sub-basin screening model. Wat 
Resour. Res. Ins. Univ. of Nebraska. 

Viessman et al, (1975). A screening model for water resources 
planning. Water Resour. Bull. 11(2), 245-255. 

Waingartner, H.K. (1967). Mathematical programming and the 

analysis of capital budgeting problems. Markham Pub. Co., 
Chicago. 

Weiss and Beard (1971). A multi basin planning strategy. 

Water Resour. Bull. 7(4), 750-764. 

Young, G. K. , Mosely, J.C. and Evenson, D, E. (1969). Sequencing 
of element construction in a multi reservoir system. 

Water Resour. Bull. 6, 528-541. 





